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1 Introduction 



Recently [1] we have developed the Nambu-Jona-Lasinio model of light nuclei [2] , or 
differently the nuclear Nambu-Jona-Lasinio (NNJL) model, invented for the description 
of low-energy nuclear forces at the quantum field theoretic level. We have shown that 
the NNJL model is fully motivated by QCD [1]. The deuteron appears in the nuclear 
phase of QCD as a neutron-proton collective excitation, a Cooper np-pair, caused by a 
phenomenological local four-nucleon interaction. Strong low-energy interactions of the 
deuteron coupled to itself and other particles are described in terms of one-nucleon loop 
exchanges. The one-nucleon loop exchanges allow to transfer nuclear flavours from an 
initial to a final nuclear state by a minimal way and to take into account contributions 
of nucleon-loop anomalies determined completely by one-nucleon loop diagrams. The 
dominance of contributions of nucleon-loop anomalies is justified in the large N c approach 
to the description of non-perturbative QCD with SU (N c ) gauge group at Nc — > oo, where 
Nc is the number of quark colours. 

In this paper we apply the NNJL model to the description of low-energy nuclear forces 
for electromagnetic and weak reactions with the deuteron by example of the evaluation of 
the cross sections for the neutron-proton radiative capture n + p — > D + 7 for thermal 
neutrons caused by the 1 S — > 3 Si transition (the Ml-transition), the photomagnetic 
disintegration of the deuteron 7 + D — > n + p and weak reactions of astrophysical interest: 
1) the solar proton burning p + p — > D + e + + z/ e , 2) the pep-process p + e~ + p — > 
D + v c and 3) the reactions of neutrino and anti-neutrino disintegration of the deuteron 
caused by charged z/ e + D — >c~+p + p, z/ c + D^e + + n + n and neutral z/ e (z/ e ) + D 
— > v c (u c ) + n + p weak currents. 

Low-energy electromagnetic nuclear reactions with the deuteron. It is well 

known that the reaction of the neutron-proton radiative capture n + p — > D + 7 for 
thermal neutrons plays an important role for the primordial nucleosynthesis in the Big 
Bang model [3]. Indeed, the deuterons produced via the neutron-proton radiative capture 
n + p — > D + 7 burn to 4 He through the reactions D + p — > 3 He + 7 and 3 He + n — > 4 He 
+ 7. Therefore, the correct description of the neutron-proton radiative capture n + p — > 
D + 7 is essential for the theoretical prediction of the amount of the matter in Universe. 

The reaction of the photomagnetic disintegration of the deuteron 7 + D — > n + p 
is related to the neutron-proton radiative capture n + p — > D + 7 via time-reversal 
invariance of strong and electromagnetic forces. The investigation of this reaction is 
conceived to get an additional check of our result for the Ml-capture related to the 
analysis of the energy dependence of the cross section calculated in the NNJL model at 
energies far from threshold. 

The cross section cr(np — > Dj)(T Q ) for the neutron-proton radiative capture has been 
measured for thermal neutrons at the laboratory kinetic energy T n = 0.0252 eV that 
corresponds to the laboratory velocity v n /c = 7.34 x 10~ 6 (the absolute value is v n = 
2.2 x 10 5 cm/s) [4]: 

<7(np -> D 7 ) cxp (T n ) = (334.2 ± 0.5) mb. (1.1) 

For the first time the cross section for the neutron-proton radiative capture has been 
calculated by Austern [5] in the Potential model approach (PMA) : 

^(np -> D 7 )(T n ) = (303 ± 4) mb. (1.2) 
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The observed discrepancy about 10% has been then explained by Riska and Brown [6] in 
terms of the contributions of the exchange currents and the A(1232) resonance. 

Recently the evaluation of the cross section for the neutron-proton radiative capture n 
+ p — > D + 7 has been carried out by using Chiral perturbation theory in the framework 
of the Effective Field Theory (EFT) approach [7-10] formulated by Weinberg within 
Effective Chiral Lagrangian description of nuclear forces [11] (see also Refs. [12]). The 
theoretical results obtained in Refs. [7-10] for the cross section for the neutron-proton 
radiative capture are rather contradictory. Indeed, in Ref. [7] the experimental value 
of the cross section for the Ml-capture has been reproduced without free parameters, 
the definition of which demands a fit of experimental data, with an accuracy better 
than 1%, a(np — > D^y)(T n ) = (334 ± 3)mb. In turn, in the more recent publications 
[8-10] the predictions are not so much optimistic. Indeed, in Ref. [8] there has been 
obtained the value a(np — > Dj)(T n ) = 297.2 mb which is about 11% less compared with 
the experimental one Eq. (|l.l|) . As has been stated in Ref. [8] for the evaluation of the 
correct value of the cross section for the Ml-capture within the EFT one needs to add 
a free parameter undefined in the approach. This parameter should be fixed from the fit 
of the experimental value Eq. Ql.lj ). This program has been realized in Ref. [9]. Then, in 
Refs. [10] the neutron-proton radiative capture has been calculated for the center of mass 
energies of the np pair up to T np < 1 MeV, T n = 2 T np , by including the contribution of the 
El-transition in addition to the Ml. This has added new free parameters with respect to 
that introduced in Refs. [8,9], where only the Ml-transition has been taken into account. 

In the NNJL model we calculate the cross section for the Ml-capture both in the 
tree-meson approximation and by including the contributions of chiral one-meson loop 
corrections and the A (1232) resonance. For the evaluation of chiral one-meson loop 
corrections we apply Chiral perturbation theory at the quark level (CHPT) 9 [13] developed 
within the extended Nambu-Jona-Lasinio (ENJL) model with a linear realization of chiral 
£7(3) x i/(3) symmetry. 

In our consideration the A(1232) resonance is the Rarita-Schwinger field [14] A" (a?), 
the isotopical index a runs over a = 1, 2, 3, having the following free Lagrangian [15,16]: 

££ n (x) = A£(x)[-(*7 Q <9 Q - M A ) sT + ^tVC*/ ^ - M A )^Y\K(^), (1-3) 

where Ma = 1232 MeV is the mass of the A resonance field A" (a;). In terms of the 
eigenstates of the electric charge operator the fields A"(x) are given by [15,16] 




The fields A* (a;) obey the subsidiary constraints: <9 M A^(x) = 7 At A^(x) = [14-16]. The 
Green function of the free A-field is determined by 

(0|T(A M (x 1 )A,(x 2 ))|0) = -iS^fa - x 2 ). (1.5) 

In the momentum representation S fJill (x) reads [15-17]: 

SM = M^Tp [ - ^ + 3 W + 3 M A + 3 Mi ) " (L6) 
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The most general form of the 7rNA- interaction compatible with the requirements of chiral 
symmetry reads [18]: 




-A+(x)e^p(x)dy(x) - A° (x)G^p(x)d ip Tr°(x) + ... , 



(1.7) 



where — M n = M p = 940 MeV is the nucleon mass. The nucleon field N(x) is the 
isotopical doublet with the components N(x) = (p(x),n(x)), and ir a (x) is the pion field 
with the components ^(x) = (n~(x) + n + (x))/\/2, n 2 (x) = (tt~(x) — n + (x))/i\/2 and 
vr 3 (x) = n°(x). The tensor is given in Ref. [15]: W = g^-(Z+l/2)^^, where the 
parameter Z is arbitrary. There is no consensus on the exact value of Z. From theoretical 
point of view Z — 1/2 is preferred [15]. Phenomenological studies give only the bound 
\Z\ < 1/2 [18]. The value of the coupling constant svna relative to the coupling constant 

9-n-NN IS g-^NA = 2(? vr NN [19]. 

Assuming that the transition A — > N + 7 is primarily a magnetic one the effective 
Lagrangian describing the A — > N + 7 decays can be determined as [19-21]: 



where F a/3 (x) = d a A l3 ( y x)—d l3 A a (x) is the electromagnetic strength field tensor and A a (x) 
is a vector potential of the electromagnetic field. The value of the coupling constant # 7 na 
relative to the coupling constant svnn is S^na = 0.14 ^nn caused by the SU (6) symmetry 
of strong low-energy interactions [19]. 

The NNJL model realizes the Lagrange approach to the description of low-energy 
nuclear forces [1]. For the evaluation of the effective Lagrangian of the transition n + p 
— > D + 7 it is necessary, first, to determine the effective Lagrangian of the strong low- 
energy transition n + p — > n + p, or more generally N + N — > N + N, where N = (p, n) is 
a nucleon field. Since the NNJL model describes low-energy interactions of the deuteron 
in terms of one-nucleon loop exchanges the effective Lagrangian of the transition N + N 
— > N + N plays an important role. Due to the transition n + p — > n + p the np pair 
on-mass shell in the initial state transfers itself into the np pair off-mass shell couples to 
the deuteron and the photon through one-nucleon loop exchanges. Then, the one-nucleon 
loop diagrams are calculated at leading order in the 1 / Mn expansion that corresponds to 
the large N c expansion due to the proportionality M N ~ AVj valid in the multi-colour 
QCD with SU(N C ) gauge group at N c -> 00 [22]. 

Such a procedure of the evaluation of effective Lagrangians in the NNJL model resem- 
bles that has being used in the ENJL model [13,23-25] for the derivation of effective chiral 
Lagrangians up to the formal replacement q(q) — ■> N(N), where q(q) is a quark (anti- 
quark) field. In the ENJL model the dominance of the leading order contributions in the 
1/Mg expansion, where M q is a constituent quark mass, has been explained by a dynamics 



£ 7 na(^) = ie 



^6 M N 



le # 7 na 



^N(xh al 5 A*(x)F^(x)+h.c.= 
\p(x) lal ^+{x) + n{x) lal ^l{x)} F? a (x) + h.c, 



(1.8) 
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of quark confinement, whereas in the NNJL model the dominance of the leading contri- 
butions in the 1/Mn expansion is justified by the large Nc approach to non-perturbative 
QCD. 

Since relative momenta of the np pair in the reaction of the Ml-capture are smaller 
compared with the mass of the pion M n = 135 MeV, for the derivation of the effective 
Lagrangian of the strong low-energy transition n + p->n + pwe can follow the ideology 
of the EFT [7-12] and integrate out pion degrees of freedom as well as other heavier 
degrees of freedom. The result of the integration can be represented in the form of two 
contributions, where the first is the explicit one-pion exchange, whereas the second is a 
phenomenological one describing a collective contribution of both many-pion exchanges 
and heavier meson degrees of freedom such as er(660), p(770), c<j(782) and so on. 

The effective Lagrangian C^^ ion (x) of the strong low-energy transition n + p — > n 
+ p in the one-pion exchange approximation is defined by 

CV-ZnW = % {[n(x) 7 5 n(x)][p(x) 7 5 p(x)] - 2 \p(x)>fn{x)][n{x)>fp(x)]}, (1.9) 

where svnn = 13.4 is the 7rNN coupling constant, n(x) and p(x) are the operators of 
the interpolating fields of the neutron and the proton. The effective Lagrangian Eq.(L9) 



is local, as we have neglected the squared momentum transfer — q 2 with respect to the 
squared pion mass, — q 2 <C M 2 . Since in the reaction n + p — > D + 7 the np pair couples 
to the deuteron and the photon in the 1 Sq state, we should rearrange the operators of the 



neutron and the proton interpolating fields in the effective interaction Eq. ( |1.9|) by such 
a way to introduce the products of the np operators creating the np states with definite 
total spins. This rearrangement can be carried out by means of the Fierz transformation 
that gives 

+ \a^C®Ca tlVl (1.10) 

o 

where a^ u = (7 »r - 7"7 At )/2- By virtue of Eq. flTTol) we recast the effective Lagrangian 
Eq.( |1.9| ) into the form 



CZ^LW = 9 f^{[P(xh 5 n c (x)][n^xh 5 p(x)] + [pfxl/Mx)]^^)] 

+ [p(x)n c (x)][n c (x)p(x)} + 3 [p(x) r y fJ 'n c (x)}[n c (x)j tl p(x)} 
+^ \p{x)a^n c {x)}[n c {x)a^p{x)]}. (1.11) 

Here n c (x) = n T (x)C and n c (x) = Cn T (x), where C is a charge conjugation matrix and 



T is a transposition. The first two terms in the effective Lagrangian Eq . (|1 -11 ) describe 



the strong low-energy 1 S — > 1 S transition of the np pair in the 1 S state. Thereby, the 
effective Lagrangian providing the 1 So — > x So transition of the np pair and caused by the 
one-pion exchange we would use in the form 

*~"one — pion \ ^ / 

= ^ {[p(x)^n c (x)}[nc{x)^p{x)} + ^n c (i)] [n^x) lfll 5 p(x)]}. (1.12) 
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The phenomenological part of the effective Lagrangian responsible for the strong low- 
energy 1 S — > 1 S transition of the np pair in the 1 S state we would choose by following 
the EFT ideology [7-12] as well and write 

£pT np (*) = 

= 77^ {[P(xh 5 n c (x)}[n c (x)-f 5 p(x)} + [p(x)^ 5 n c (x)][n c (x)^ 5 p(x)}}, (1.13) 

where a np = (—23.75 ±0.01) fm is the S-wave scattering length of the elastic np scattering 
in the 1 S state [17]. 

The appearance of the S-wave scattering length for the definition of the phenomeno- 
logical coupling constant of the strong low-energy transition n + p — ► n + p, or differently 
low-energy elastic np scattering is rather natural in the EFT. Indeed, in the EFT pions 
are treated perturbatively and at leading order in the pion-exchange approximation for 
the description of the S-wave scattering length of low-energy elastic np scattering in the 
1 So state Weinberg has introduced a phenomenological four-nucleon interaction with a 
coupling constant [11] 

C„ = -^f. (1.14) 

As has been stated in the EFT this constant is a collective contribution coming from the 
integration over all meson degrees of freedom heavier than the pionic ones. Since in our 
approach we separate the pionic degrees of freedom from that with masses heavier than the 
pion as well, the integration over these heavy meson degrees of freedom should have the 
same form as it is postulated in the EFT. We have only halved the coupling constant Co 
defined by Eq. ( |1.14| ) in order to distribute symmetrically phenomenological contributions 
between couplings [p(x)7 5 n c (a;)][n c (x)7 5 p(x)] and \p(x) r y flr y 5 n c (x)}[n c (x)'y fl r y 5 p(x)]. 

The total effective Lagrangian describing the strong low-energy transition n + p — > n 
+ p of the np pair in the 1 S state is then determined by 

*"eff K^J *-one-pionV J '/ ~ '"ph. \^ J 

= CW {[p(x)rfn c (x)][n c (x)rfp(x)} + \p(x)^ 5 n°(x)] [n c (x)'y fi y 5 p(x)]}, (1.15) 

where the effective coupling constant Cnn of the strong low-energy transition n + p — ► n 
+ p is equal to 

Cnn = lb - ICT* = 3 ' 27 x 10 " 3MeV " 2 - (1-16) 

Note that the contribution of the phenomenological part to the effective coupling constant 
Cnn makes up less than 33%. 

Since nuclear forces are isotopically invariant [26], the effective Lagrangian of the 
strong low-energy N + N — > N + N transition of the NN pair in the 1 So state can be 
defined as follows 

{[p(x)^ b n c (x)][n c (x)^p(x)} + [p(x)Yl 5 n c (x)][n c (x)'y^ 5 p(x)] 
+ -[n(x) , y 5 n c (x)][n c (x)'y 5 n(x)] + [n(x)^^rfn c (x)][n c (x)^^n(x)} 

+ \{p{x)^p c {x)}{p c {x)^p{x)} + [p(x)^ 5 p c (x)][p c (x)^ 5 p(x)]}. (1.17) 
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We would like to emphasize the effective Lagrangian Eq.( 1.16Q describing strong low 



energy N + N — ► N + N transitions is obtained in complete agreement with the EFT 
ideology. 

In the low-energy limit the effective local four-nucleon interaction Eq. ( |1.17| ) vanishes 
due to the reduction 

[N(x)^^ 5 N c (x)}[N c (x)^N(x)} -> -[N(x)^N c (x)}[N c (x)^N(x)}, (1.18) 

where N(x) is an operator of the neutron or the proton interpolating field. Such a vanish- 
ing of the one-pion exchange contribution to the NN potential is well-known in the EFT 
approach [11,12] and the PMA [26]. In power counting [11,12] the interaction induced 
by the one-pion exchange is of order 0(p 2 ), where p is a relative momentum of the NN 
system. The former is due to the Dirac matrix 7 5 which leads to the interaction between 
small components of Dirac bispinors of nucleon wave functions. 

In the one-nucleon loop exchange approach the contributions of the interactions 
[iV(a;)7 A1 7 5 A^ c (x)][A^ c (x)7^7 5 A^(a;)] and [A(a;)7 5 N c (x)] [A c (x)7 5 iV(:r)] to the amplitudes 
of nuclear reactions are different and do not cancel each other in the low-energy limit 
due to the dominance of nucleon-loop anomalies [1]. This provides the interaction be- 
tween large components of Dirac bispinors of nucleon wave functions that distinguishes 
the NNJL model from the EFT. 

Low— energy weak nuclear reactions with the deuteron. The weak nuclear reaction 
p + p^D + e + + z/ e , the solar proton burning, plays an important role in Astrophysics 
[3,27]. It gives start for the p-p chain of nucleosynthesis in the Sun and main-sequence 
stars [3,27]. In the Standard Solar Model (SSM) [28] the total (or bolometric) luminosity 
of the Sun L Q = (3.846 ± 0.008) x 10 26 W is normalized to the astrophysical factor 5 , pp (0) 
for the solar proton burning. The recommended value 5 , pp (0) = 4.00 x 10~ 25 MeVb [29] has 
been found by averaging over the results obtained in the Potential model approach (PMA) 
[30,31] and the Effective Field Theory (EFT) approach [32,33]. As has been shown recently 
in Ref. [34] the inverse and forward helioseismic approach confirm the recommended value 
of 5^(0) within experimental errors on the helioseismic data and solar neutrino fluxes. 

In this paper we apply the NNJL model to the description of low-energy nuclear forces 
for weak nuclear reactions with the deuteron of astrophysical interest: 1) the solar proton 
burning p + p — > D + e + + z/ e , 2) the pep-process p + e~+p— >D + z/ e and 3) the 
reactions of neutrino and anti-neutrino disintegration of the deuteron caused by charged 
f e + D— > e~ + p + p, u c + D-^e + + n + n and neutral v e {p e ) + D — > u c (u c ) + n + p 
weak currents. The reactions z/ e + D — >e - +p + p and is e + D — > i/ e + n + p caused 
by charged and neutral weak currents, respectively, and induced by solar neutrinos are 
planned to be measured for solar neutrino experiments at Sudbury Neutrino Observatory 
(SNO) [35]. In turn, the cross sections for the reactions u c + D—>-e + +n + n and z/ e + D 
— > u c + n + p caused by charged and neutral weak currents, respectively, and induced by 
reactor anti-neutrinos have been recently measured by the Reines's experimental group 
[36]. In the sense of charge independence of weak interaction strength the observation of 
the reaction v e + D — > e + + n + n is equivalent to the observation of the reaction of the 
solar proton burning p + p — > D + e + + u e in the terrestrial laboratories. 

The paper is organized as follows. In Section 2 we evaluate the amplitude of the 
Ml-capture n + p^D + 7in the tree-meson approximation. The contribution of 
low-energy elastic np scattering to the amplitude of the process n + p— ► D + 7 is 
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obtained in agreement with low-energy nuclear phenomenology. In Section 3 we evaluate 
contributions of chiral one-meson loop corrections to the amplitude of the Ml-capture 
in Chiral perturbation theory at the quark level (CHPT) g developed within the ENJL 
model with a linear realization of chiral U(3) x U (3) symmetry. In Section 4 we include the 
contribution of the A(1232) resonance and analyse the total cross section for the neutron- 
proton radiative capture for thermal neutrons and compare it with experimental data. In 
Section 5 we treat the reaction of the photomagnetic disintegration of the deuteron 7 + 
D — > n + p related to the neutron-proton radiative capture n + p — > D + 7 via time- 
reversal invariance and analyse the energy dependence of the cross section at energies far 
from threshold. In Sect. 6 we evaluate the amplitude of the solar proton burning. We 
show that the contribution of low-energy elastic pp scattering in the 1 So state with the 
Coulomb repulsion is described in agreement with low-energy nuclear phenomenology in 
terms of the S-wave scattering length and the effective range. In Sect. 7 we evaluate the 
astrophysical factor for the solar proton burning and obtain the value S pp (0) = 4.08 x 
10~ 25 MeVb agreeing good with the recommended one S pp (0) = 4.00 x 10 _25 MeVb. In 
Sect. 8 we evaluate the cross section for the neutrino disintegration of the deuteron v c 
+ D^e~+p + p with respect to S , pp (0). In Sect. 9 we adduce the evaluation of the 
astrophysical factor S , pep (0) for the pep-process relative to S , pp (0). In Sects. 10 and 11 we 
evaluate the cross sections for the anti-neutrino disintegration of the deuteron P e + D — > 
e + + n + n and z/ c + D^z/ c + n + p and average them over the anti-neutrino energy 
spectrum. The average values of the cross sections agree well with experimental data [36]. 
The cross sections for the weak nuclear reactions of astrophysical interest are calculated 
at zero contribution of the nucleon tensor current [1]. This makes the description of 
low-energy nuclear forces within the NNJL model compatible with the predictions of 
the SSM [27,28]. In more detail this is discussed the Conclusion and Appendix. In 
the Conclusion we discuss the obtained results. In Appendix we evaluate the effective 
Lagrangian C^~^ Dc Ue (x) of the low-energy weak transition p + p — > D + e + + u e . The 
contribution of the nucleon tensor current [1] to the effective Lagrangians of low-energy 
weak nuclear transitions of astrophysical interest like p + p^D + e + + z/ e and so on is 
evaluated and discussed. 



2 The Ml— capture in the tree— meson approximation 

Since the NNJL model realizes a Lagrange approach to the description of low-energy 
nuclear forces [1], first thing what we have to do is to evaluate the effective Lagrangian 
£°fj^ D7 (x) of the transition n + p — > D + 7. In the tree-meson approximation the effective 
Lagrangian £™^ Dl (x) is defined by one-nucleon loop diagrams depicted in Figs. 1 and 2. 
The evaluation of these diagrams at leading order in the large N c expansion yields 

+ ^p-/O^0C NN M N I}1^ (2.1) 

where *F^ u (x) = \ e^ af3 F a p(x), /x p = 2.793 and /i n = —1.913 are the magnetic dipole 
moments of the proton and the neutron measured in nuclear magnetons, D^(x) is the 
operator of the interpolating field of the deuteron and D^ v (x) = d fl D u (x) — d^D^x), 
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5V is a phenomenological coupling constant of the NNJL model related to the electric 
quadrupole moment of the deuteron Q D = 0.286 fm 2 [1]: g\ = 2tc 2 Q d M^. 

The matrix element of the transition n + p^D + 7we define by a usual way 

d 4 x<z>(fc D hWI4r D >)MPi)p(p 2 )> = 
W^tb+i-a-a) + + (22) 



where i?j (i = 1, 2, D) and are the energies of the neutron, the proton, the deuteron and 
the photon, V is the normalization volume. 

The wave functions of the initial \n(pi)p(p 2 )) and final (D(k-[))^f(k)\ state we take in 
the usual form 

\n(pi)p(p 2 )} = a{(p 1 ,a 1 )al(p 2 ,cr 2 )\0}, 

(D{k D )j{k)\ = (0\a D (k D ,\ D )a(k,\), (2.3) 

where a^(pi, &i) and aUp 2 , cr 2 ) are the operators of creation of the neutron and the proton, 

and an(kT), Ad) and a(k,X) are the operators of annihilation of the deuteron and the 
photon. 

The matrix element of the transition n + p — > D + 7 reads 

A*(n + p-D + 7 ) = (^ p -^ n )^-^C' NN e^ a e;(A ; ,A)e^D,AD) 

x[u c (p 2 )(2k Du - M N7 ,)7 5 M ( Pl )], (2.4) 

where e%(k, A) and e* (kn, Ad) are the 4-vectors of the polarization of the photon and the 
deuteron, then u c (p 2 ) and u(pi) are the bispinorial wave functions of the proton and the 
neutron, respectively, normalized by u c (p 2 )u c (p 2 ) = — 2M N and u(pi)u(pi) = 2M N . 
In the low-energy limit when 

[u c {P2hulMPi)} -> ~9uo [u c {p 2 h 5 u{Pi)] (2-5) 
and k-D U —>■ g u o2M-^ the matrix element Eq. (|2.4j ) acquires the form 



A<(n + p^D + 7) = e(/Xp-/x n )^C , N N(^Xe*(^,A))-e*(^D,A D ) 

x[« c (p 2 )7MPi)]- (2.6) 

The evaluation of the matrix element Eq. (|2.6|) , the effective vertex of the n + p — > D + 
7 transition, we have carried out with the wave functions of the neutron and the proton 
in the form of the plane waves. However, a physical 1 So state of the np pair is defined by 
low-energy nuclear forces. For the description of the contribution of low-energy nuclear 
forces to the physical 1 S state of the np pair coupled to the deuteron and the photon 
we suggest to sum an infinite series of one-nucleon bubbles with vertices defined by the 
effective Lagrangian £°g~* np (x) Eq.( |1.15| ). The result of the summation can be represented 
in the following form 

M(n + p -> D + 7) = e (jip - fi n ) C NN (k x e*(k, A)) ■ e*(k D , A D ) [u c (p 2 )-f 5 u( Pl )] 

16tt 2 J ixH \ M N -q-P-Q M N - q - Q J 
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where P = pi + P2 = (2yp 2 + Mj|, ) is the 4-momentum of the np pair in the center 
of mass frame. Then, Q = a P + b K = a (p± + P2) + b (pi — P2) is an arbitrary shift of 
virtual momentum with arbitrary parameters a and b, and in the center of mass frame 
K = pi — p 2 = (0,2p). The explicit dependence of the momentum integral on Q can be 
evaluated by means of the Gertsein-Jackiw procedure [37] (see also Ref. [1]). It is given 
by 

A tr f . 1 U 

7rh { Mn — q — P — Q M n — q — Q) 

tr{ 7 5 . 7 5 — —\ -2a(a + l)P 2 -2b 2 K 2 . (2.8) 

vr 2 « \ M N - q - P M N - q J V ; V 1 

For the evaluation of the momentum integral over q we would keep only the leading order 
contributions in the 1/Mn expansion caused by the large Nc expansion [1]. This yields 

A f 5 1 5 1 1 

-^tr 7 5 - 3-^7 5 1 



vr 2 2 {' M N -q-P-Q' M N — q — Q . 
-8 a (a + 1) M N + 8 (6 2 - a (a + 1)) p 2 — i8%M N p. (2.9) 



The expression Eq. (|2.7| ) we reduce to the form 
.M(n + p^D + 7 ) = e(/i p -/i n )^C , NN(A?xe*(fc,A))-e*(fc D ,A D )[n c (p 2 )7 5 M(p 1 )] 



2 



1 2 

1 - -r Q pa n pP + ia np p 



(2.10) 



Here we have denoted 



C NN M N 2 2 1 

fl n P = - Z , r np = (b -a(a+l))- — , 

2ir 7r Mn 

1 = !_^(|±i)c NN MS, (2.11) 

where r np = 2.75 ± 0.05 fm is the effective range of low-energy elastic np scattering [17]. 

Renormalizing the wave functions of nucleons \/Zu{px) — > u(pi) and \/~Zu(p2) — > 
■u(p 2 ) we arrive at the expression 

.M(n + p^D + 7 ) = e(/i p -/i n )^C N N(fcxe*(fc,A))-e*(fc D ,A D )[u c (p 2 )7 5 M(pi)] 

x j , (2.12) 

1 - -r np a np p 2 + i a np p 

where the factor 1/(1 — \r ap a np p 2 +i a np p) describes the contribution of low-energy nuclear 
forces to the physical 1 So state of the np pair coupled to the deuteron and the photon. 
It has the form of the amplitude of low-energy elastic np scattering in the 1 S state in 
complete agreement with low-energy nuclear phenomenology [26]. By using the relation 
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expressing the phase shift 5 np (p) of low-energy elastic np scattering in terms of the S-wave 
scattering length a np and the effective range r np [26] 

ctg 5 np (p) = - — + l-r np p 2 (2.13) 
we can recast the factor 1/(1 — ^r np a np p 2 + ia np p) into the form 



1 ~ ^r np a np p 2 + i a np p anp 

In terms of the phase shift 5 np (p) the expression Eq. fl2.12| ) reads 

A-l(n + p ->D + 7) = e (// p - /i n ) CW x e*(k, A)) -e*(k D ,\ D ) 

x [uc(p 2 h 5 u(Pi)} e^p(P) Sin * np(p) . (2.15) 

-a np p 

In the NNJL model low-energy nuclear forces between the neutron and the proton in 
the physical deuteron state are described by the one-nucleon loop exchanges in terms of 
the phenomenological coupling constant <jv which is defined by the electric quadrupole 
moment of the deuteron Q D , g\ = 2tt 2 Q d M^ [1]. The electric quadrupole moment of the 
deuteron is caused by nuclear tensor forces [26]. Therefore, the relation g\ = 2tt 2 Q d M^ 
confirms at the quantum field theoretic level the fact pointed out by Blatt and Weisskopf 
that the existence of a bound triplet state of the neutron-proton system would be entirely 
due to the tensor force [38]. Thus, in NNJL model through the phenomenological coupling 
constant gy tensor forces govern the existence of the deuteron as a bound neutron-proton 
triplet spin state and strength of low-energy interactions of the deuteron with nucle- 
ons and other particles in terms of the one-nucleon loop exchanges. The evaluation of 
one-nucleon loop diagrams at leading order in the 1/Mn expansion, or in the large Nq 
expansion, reduces a momentum dependence of nucleon diagrams to the trivial form ac- 
counting for only the Lorentz covariant properties of the interaction. In this approach the 
deuteron looks like a point-like particle. Such a representation is enough for the evalua- 
tion of effective Lagrangians of different low-energy nuclear transitions with the deuteron 
in an initial or a final state describing effective vertices of low-energy nuclear transitions 
defined at their thresholds. However, for the evaluation of amplitudes of low-energy nu- 
clear reactions for energies far from threshold one needs to take into account a spatial 
smearing of the physical deuteron caused by a finite radius r D = 1/V^d^n = 4.319 MeV 
[17] determined by the binding energy of the deuteron 6d = 2.225 MeV. The spatial 
smearing of the physical deuteron we introduce phenomenologically in the form 

MP) = ; \ 2 , (2-16) 
1 + r^p z 

that is nothing more than the momentum representation of the approximate 3 Si wave 
state of the deuteron [26] . 



= e z5 np (p)^M. (2.14) 
-a np p 
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Substituting Eq. (|2.16|) in Eq. (|2.15|) we obtain the amplitude of the Ml-capture calcu- 
lated in the NNJL model: 



.M(n + p^D+ 7 ) = e(p p - p n )-^C NN (kx e*(k,\))-e*(k D ,\ D ) 

x [u^hMpi)} A(P)^!-L^. (2.17) 

-a np p 1 + rfyp 2 

For thermal neutrons the kinetic energy of the relative movement of the np pair is of 
order T np ~ 10~ 8 MeV [9]. This yields the relative momentum of the np pair to be smaller 
compared with the binding energy of the deuteron, p~3x 10 -3 MeV <C £d = 2.225 MeV. 
Therefore, for thermal neutrons without loss of generality we can calculate the amplitude 
of the Ml-capture setting p = 0: 

_M(n + p^D + 7 ) = e(/i p -/i n )^C NN (fcxe'(M)).e , (fc D ,A D ) 

x [uc(p 2 )^u( Pl )]. (2.18) 

Thus, we have obtained that the amplitude of the neutron-proton radiative capture n 
+ p — > D + 7 for thermal neutrons coincides with the matrix element of the effective 
Lagrangian C^~^ D,y (x) given by Eq.( |2.6| ). 

The cross section for the Ml-capture for thermal neutrons calculated in the tree-meson 
approximation is then defined by 

a(np^ D 7 )(T n ) = 1 _ ^ n ) 2 g ^Q D C 2 NN M N + 1 ^ = 276mb. (2.19) 

The theoretical value cx(np — > D7)(T n ) = 276 mb is about 17% smaller compared with 
the experimental one <r(np — > D7) exp (T n ) = (334.2 ± 0.5) mb. Thus, in the tree-meson 
approximation the NNJL model predicts the cross section for the Ml-capture for thermal 
neutrons somewhat worse than the EFT, cx(np — > D7)(T n ) = 297.2 mb [8]. In order to 
improve the agreement with the experimental data we have to include chiral one-meson 
loop corrections [7-9] and the contribution of the A(1232) resonance [6,7]. 

3 Chiral one— meson loop corrections to the ampli- 
tude of the Ml— capture 

For the evaluation of chiral meson-loop corrections in the NNJL we use (CHPT) g 
developed in Refs. [13] within the ENJL model with a linear realization of chiral U(3) x U(3) 
symmetry. Below we consider the contributions of chiral one-meson loop corrections 
induced by the virtual meson transitions n — > a^, a\ — > ttj, n — > (cu,p)7, (uj,p) — > 
7T7, a — > (u, p)j and (ui, p) — > aj, where a is a scalar partner of pions under chiral 
SU(2) x SU (2) transformations in (CHPT)q with a linear realization of chiral U (3) x £7(3) 
symmetry [13]. 

The effective Lagrangians 5C^{x) and 5C^P(x), caused by the virtual meson tran- 
sitions 7r — > ai 7, a\ — > 7r 7, 7r — > (ou,p)j, (u>,p) — > ttj , a — > (u , p)j and (u , p) — > aj, 
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we evaluate at leading order in the large Nq expansion [1]. The results of the evalua- 
tion contain divergent contributions. In order to remove these divergences we apply the 
renormalization procedure developed in (CHPT) 9 for the evaluation of chiral meson-loop 
corrections (see Ivanov in Refs. [13]). Since the renormalized expressions should vanish 
in the chiral limit M n — > [13], only the virtual meson transitions with intermediate 
7r-meson give non-trivial contributions. The contributions of the virtual meson transi- 
tions with the intermediate cx-meson are found finite in the chiral limit and subtracted 
according to the renormalization procedure [13]. Such a cancellation of the cr-meson con- 
tributions in the one-meson loop approximation agrees with Chiral perturbation theory 
using a non-linear realization of chiral symmetry, where a-meson-like exchanges can ap- 
pear only in two-meson loop approximation. Then, the sum of the contributions of the 



virtual meson transitions tt 



p 7, 7T 



o 



p°7 and 7r° — > to the effective coupling mry 



is equal to zero. As a result the effective Lagrangians 5C^p(x) and S£ 



nri7 / 



are given by 



5C 



PP7, 
eff 1 



X 



a r 



9a9ttNN 



16tt 3 

x \p(x)a^p(x)} F^(x) 



M n J^N + S'ttNN 



N c a n M\ 



' CU ' MN Ml J. 



16vr 3 F w 



ttVN 



ie 
4M 



— 9a9ttNN 



at. 



16tt 3 F, 



J; 



7raiN 



n(x)a pv n(x)]F^(x), (3.1) 



where a p = g 2 „/^ = 2.91 is the effective coupling constant of the p — > decay, F n 



92.4 MeV is the leptonic coupling constant of pions, and g& 
andJ^vN are the momentum integrals determined by 

d 4 p 



1.267 [39]. Then, J, 



7raiN 



J; 



7raiN 



1 



J - 



VN 



TT 2 {Ml + p 2 ) (Ml + p 2 )(M 2 + f) 

d 4 p 1 

~^ (M 2 + p 2 ) (My + p 2 ) (M 2 + p 2 ) 



0.017M" 



0.024 M: 



M r , .\/_. 770 MeV [39] and M a 



(3.2) 
V2M P 



where p is Euclidean 4-momentum, My 

At Nc = 3 the cross section for the Ml-capture accounting for the contribution of the 
effective interaction Eq.( [3.1|) amounts to 

/ \ /m \ / \ *5 25 (y. 

cr(np -> D7)(T n ) = - (p p - p n ) 



x 



1 + 



PttNN 




' 2 OA ~ ^ D ^NN 
04 7T Z 



p p - p n 8tt 2 tt \ v " 7rai " ' 2g A 

where we have used the relation svnn — gxM^/F^. The theoretical value of the cross 
section for the neutron-proton radiative capture given by Eq. ( |3.3|) differs from the exper- 
imental one by about 14%. This discrepancy we should describe by taking into account 
the contribution of the A(1232) resonance. 



N + 



VN 



287.2 mb, 



(3.3) 



4 The A(1232) resonance 

For the evaluation of the contribution of the A(1232) resonance to the amplitude of 
the Ml-capture in the NNJL model we have to obtain the effective Lagrangian £g|^ AN (x) 
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describing the strong low-energy n + p — > A + N transition. For this aim we should use 
the procedure having been applied to the evaluation of the effective Lagrangian C^~^ np (x) 
given by Eq. fll .1 5|) . In Refs.[6,7] the evaluation of the contribution of the A(1232) res- 
onance in terms of exchange currents has been carried out in the one-pion exchange 
approximation. Thereby, following Refs. [6,7] we suppose to evaluate the effective La- 
grangian £gg~* AN (x) of the strong low-energy n + p — > A + N transition in the one-pion 
exchange approximation. This gives 

CT AN (*) = ff| / d*z ±^\z - x) {[A+(s) 0%n c (x)] 

x [n c (z)>y 5 p(x) + n c (x)>y 5 p(z)} - [A°(z) 0%p c (x)] [n c (z)>y 5 p(x) + n c (x)<y 5 p(z)} 
+1®7 5 ^-7^7V}, (4-1) 

where we have kept only the terms contributing to the transition of the np pair in the 1 So 
state into the AN state. Using then the phenomenological Lagrangian 

Cn P D(x) = -igv[p c (x)^n(x) - n c (x)^p(x)]Dl(x) (4.2) 

the effective Lagrangian describing the contribution of the A(1232) resonance to the low- 
energy transition n + p^D + 7is defined by 

J d 4 x C^ AN ^ D7 (x) =- / d A x 1 d 4 x 2 d 4 x 3 (T(£^ AN (xi)C npD (x 2 )C jNA (x 3 ))) = 

1 69v ^ NA ^ NA ^ NN ! d' Xl d 4 x 2 d% I d 4 z -^-5 {A \z - x x ) 



6 Mn SVnn 9-kNN 4M2 J J dz v 

x T([p c (x 1 ) 7 5 n(2) +p-(z)^n(x 1 )] D\ (x 2 )F a(3 (x 3 )) 

x { (0|T([A+(*) Q% n c ( Xl )] [p c (x 2 )Yn(x 2 ) - n c (x 2 ) 7 >(x 2 )] [p(x 3 ) lf5l 5 At(x 3 )]) |0) 
-(0|T([A° (z) Q%p c (x)][p^x 2 )Yn(x 2 ) - n^x 2 )^p(x 2 )] 
x [n{x 3 ) lpl 5 A° a (x 3 )})\0) + ( 7 5 ® 1 - - 7 ,7 5 ® 7")} = 

i e 9v S'-yNA A^nn f ,a ,a ,a f ,a ® sU)r \ 

= q 77 — 7777 / d Xld X2d x 3 d z >(z - x{) 

3 M N svnn AM* J J dz v 

x T([j9 c (x 1 ) 7 5 n( 2 ) +p c (^) 7 5 n(x 1 )] J Dt( X2 ) F ^( X3 )) 

x {(0|T([A+(z) 0%n c (x 1 )][^(x 2 ) 7 >(x 2 )][p(x3)7^7 5 A+(x 3 )])|O) 

+ (O|T([AO(^)0%p c (x 1 )]^(x 2 ) 7 ^(x 2 )][n(x 3 )WA°^3)])|O) 

+ ( 7 5 ® 1 -> -^ul b ®1 U )} = 

= |^^^f%| / d*x^x 2 d*x 3 ( d^z ±iP\z - Xl ) 
3 g wNN svnn 4M^ 7 7 9^ 

x {T([^(x 1 ) 7 5 n(z) + p fc (z) 7 6 n(zi)] £)t (x 2 )F^(x 3 )) 
x T tr{5 , Q£J (x 3 - z) 6% - x 2 )j**S F (x 2 - x 3 ) 7/37 5 } 

-T([pfe(ar 1 )7 I/ 7 5 n(z) + p c (z) lvl 5 n(x 1 )} £>t {x 2 )F^{x 3 )) 

x -tr{^(a; 3 - z) W^-fS^ - x 2 )-fS F (x 2 - x 3 ) 7/37 5 }}. (4.3) 
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Thus, the effective Lagrangian £"^ AN ^ D7 (a;) is equal to 

f j4 / -np-+AN^D7/ \ 2 ieg\ g-xNA g^yNA 9wNN f iA t4 j4 f jA ® r(4) / \ 

J dxCI {x) = 3 Mf^^iMf y d X2rf X3 i )( ^ 

x {T([p c (x 1 ) 7 5 n( 2 ) +p c (z) 1 5 n(x 1 )]Dl(x 2 )F^(x 3 )) 
x -tr-fS^a^ - z) 6% Sp(x 1 - x 2 )YS F {x2 - x 3 )^ 5 } 
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-T([p c (x 1 ) 7i , 7 5 n(z) +p c (z) 7! , 7 5 n(x 1 )] £>t (x 2 )F a/3 (x 3 )) 

x -tr{^(x 3 - z) 6^7^(07 - x 2 )^S F (x 2 - x 3 )7/37 5 }}- (4.4) 

In the momentum representation of the baryon Green functions the effective Lagrangian 
Eq. (pL"4D reads 

,4 /,np-^AN-*D7/ \ 2 iegy g n NA 9*yNA S'irNN /" j4 /" j4 ^ r(4) / \ 

^4ff W = o!72 ? — 7777 / rf x i / d z ( z ~ x i) 



3 AMI J J dz> 



d*x 2 d 4 k 2 d*x 3 d 4 k 3 ik<2 . ( X2 _ Xl ) ik 3 . (x 3 - z) 
(2vr) 4 (2vr) 4 

x{T([p c (xi) 7 5 n(z) +p c (z) 1 5 n(x 1 )}Dl(x 2 )F aP (x 3 )) 

x f^l e ikt ■ (X! - z) tr{5 - (A . + fc ) Q 1 , ^VY} 

i 1 1 Mn — ki + k 2 Mn — ki 

-T([p c (a; 1 ) 7l , 7 5 n(z)+p c (z) 7 , 7 5 n(a; 1 )] J D; t (x 2 )F Q/3 (x3)) 

x /^ie^-^-^tr^^ + ^e^ 1 ^ VY}}- (4.5) 

•/ ?r 2 z M N - h + k 2 M N - ki J 

The effective Lagrangian Eq.( [4.5| ) defines the contribution of the A(1232) resonance to 
the low-energy transition n + p — > D + 7. 

The matrix element of the neutron-proton radiative capture caused by the contribution 
of the A(1232) resonance exchange is equal to 

M(n + p -> AN -> D + 7) = 

" * 6 9V 97rNA 9lNA 9lNN ^tVpO] (k a e* p (k) - k^ a {k)) e;(k D ) 



2M2 6tt 2 ^ nn ^ nn 4M2 
xjf a (k»,k) 



2M§ bir z 4M^ 
x J^ a (k D ,k), (4.6) 

where the structure functions J^ a {kv>-, k) and J^ a {kji, k) are defined by the momentum 
integrals 

jf a (k D ,k) = 

= I tr {(^i + feDj^Cfci + fe) - — i — y^^yy}, 

•/ 7r*« M N - A;i + A; D M N - fci 
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= / tr{(A;i + ^)" 5aw (*i + k ) e ^ ^- — ^tV}- (4-7) 

J 7T « M N - fei + A; D M N - fci 

At leading order in the large Nc expansion the structure functions Eq. (|4.7| ) read 

jf a (k D ,k) = i(z - ^Ji M N e^ aX k DX , 

jr Pa (k B , k) = Uz-l)iM% e^ au . (4.8) 



We have neglected the mass difference between the masses of the A(1232) resonance and 
the nucleon. The matrix element of the low-energy transition n + p — ► D + 7 caused by 
the A(1232) resonance contribution is equal to 



A^ A (n + p -> D + 7) 



e 9v 
2M N 4tt 2 



1 _ g \ 8 ^ttNA ff 7 NA P^NN 

,2 / 9 gv NN ^nn 4M2 



x £ a ^k a e}(k, X)e;(k B , \ D )[u c (p 2 )(2k Du - M n7 ,) 7 5 M ( Pi )] 



1 _ ^ \ 8 ff^NA ff 7 NA ff^NN 
s 2 J 9 ^nn ^nn 4M2 



(A;xe*(A;,A))-e*(fcD,Ai 



x[n c (p 2 )7 5 «(Pi)]- 
In turn, the contribution of the nucleon tensor current [1] 



5C npD (x) 



9t 
2M^ 



[p c (x)a flu n(x) - n c (x)a flu p(x)} D ] (x 



does not depend on the parameter Z and reads 



5M A (n + p -> D + 7) = e 



1 8 svna 9 7 na 



NN 



5 SV 9 (?7rNN S'ttNN 4M2 



x(fc x e*(k, A)) • e*(fc D , A D ) [u c (p 2 )-f 5 u( Pl )}. 
The coupling constants #t and gy are connected by the relation [1] 



(4.9) 



(4.10) 



(4.11) 



9t = \ -g Y + 0(l/ JN c ). 



(4.12) 



The total amplitude of the neutron-proton radiative capture for thermal neutrons reads 

8tt 2 
3 

2g~A 



M(n + p -»• D + 7) = e (// p - /i n )g^CW(& x e*(£, A)) • e*(£ D , A D ) [u c (pa)7 S «(Pi)] 

fi^NN ^ tt (i| , , ^ T I , 1 1/3 18 PttNA P 7 NA fi^NN 



X 



1 H 77-^ — I ^TraiN + 71 ^ttVN 



7T 



1 - 2Z 1 4 ^na 3 7 NA g 



ttNN 



4M2J 



/Up - /U n 5 V 8 C N n 9 ^nn fi'vrNN 4M2 

(4.13) 
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The total cross section for the neutron-proton radiative capture is then defined by 

1/ _ ^2^ « 
yflp /i n J „ 

V 64 7T Z 



<r(np -> D7)(p) = - (/i p - /i n ) 2 — — Q D C*nn m n 4 



x 




TraiN + ~ JwVN 




8 flVNA ff 7 NA P^NN 

/ip - /i n 87T 2 7T \ v " 7rtlliN ' 2# A ""'"j ' /ip - /i n 5 V 8 C N N 9 ^NN fi^NN 4M 2 



1 - 2Z 1 4 g wNA g lNA g 



3 

ttNN 



1 2 



(4.14) 



/ip — /i n C*NN 9 (/,rNN A'ttNN 4M 2 _ 

The numerical value of the cross section amounts to 

a(np^D 7 )(T n ) = 325.5 (1 + 0.246(1 -2Z)) 2 mb. (4.15) 

Thus, the discrepancy of the theoretical cross section and the experimental value Eq. ( |1.1|) 
can by described by the contribution of the A (1232) resonance. In order to fit the ex- 
perimental value of the cross section we should set Z = 0.473. This agrees with the 
experimental bound \Z\ < 1/2 [18]. At Z = 1/2 that is favoured theoretically [15] we get 
the cross section a(np — > Dj)(T n ) = 325.5 mb agreeing with the experimental value with 
accuracy better than 3%. 



5 The photomagnetic disintegration of the deuteron 

The amplitude of the photomagnetic disintegration of the deuteron 7 + D — > n + p 
is related to the amplitude of the neutron-proton radiative capture n + p — > D + 7 due 
to time-reversal invariance and reads 



5g 



M(*y + D -> n + p) = e (/i p - /i n ) ^ C NN (k x e(k, A)) • e(k D , A D ) [u(jp 2 )^u c {p x )\ 



5„,c/ 



57T 



X 



1 H 7T"^7 I <i7raiN + ^ -'ttVN 



1 1/318 (/ttNA </ 7 NA 9tNN 



/ip - /i n 87T 2 7T 



2& 



/^p - 5 V 8 C N N 9 ^nn 4M 2 



1 — 2Z 1 4 (/ttna 37NA P^nn] z 



5 np (p) Sil1 



/i P - /in CW 9 S'tNN ^nn 4M 2 J -a np p 1 + r 2 ^ 2 ' 

The cross section defined by the amplitude Eq. QS.ip is then given by 



<r( 7 D -> np)(^) = (T 



1 2 ^np^npP 



2) 2 + a 2 pp2 (1 + rgp 



2 „2\2 ' 



(5.1) 



(5.2) 



where p = y M^(uj — e^) is the relative momentum of the np pair in the 1 So state, u is 
the energy of the photon, and cr is equal to 



2 2haQv 2 3/ 2 5/2 

O"0 = (/ip - Mn) - — „ Cmnt £n M*' 



X 



1 + 



a 1 

i/vrNN 



1927T 

M 2 a p 



— ^ NN c D m N 



/ip - /i n 87T 2 7T 



2<7A 



J, 



ttVN 



1 1/318 (/,rNA <? 7 na P^nn 

/ip - /in 5 V 8 C N N 9 ^NN S^NN 4M 2 



1 - 2Z 1 4 g^A # 7 NA ^ 



NN 



/ip — /in C*NN 9 (/ttNN fi^NN 4M 2 



7.10 mb. 
17 



(5.3) 



The energy region of the dominance of the photomagnetic disintegration of the deuteron is 



ct(7D - 


-> np)(a> 


ct( 7 D - 


-> np)(a; 


<t(jD - 


-> np)(a> 



restricted by the constraint r^p — y (a; — £d)/£d <C 1 or differently > £d = 2.225 MeV. 

The numerical values of the cross section for the photomagnetic disintegration of the 
deuteron at energies to < 2 Ed = 4.45 MeV read 

= 0.358 (0.380) mb, 

aj=2.62MeV V ' 

= 0.302 (0.327) mb, 

= 0.094 (0.128) mb, (5.4) 

OJ=4.45MeV v ' ' v ' 

where in parentheses we have adduced the theoretical values obtained by Chen and Savage 
in the EFT [10]. One can see a reasonable agreement between the results obtained in the 
NNJL model and the EFT. Thus, the spatial smearing of the physical deuteron caused 
by the effective radius td and introduced in the NNJL model phenomenologically in the 
form of the wave function Eq.( |2.16D describes well the energy dependence of the cross 
section for the photomagnetic disintegration of the deuteron at photon energies far from 
threshold. Note that in the critical region of the photon energies uo < 2 £rj = 4.45 MeV 
the cross section for the photomagnetic disintegration of the deuteron calculated in the 
NNJL model falls steeper with uo than in the EFT. However, in this energy region the 
dominant role is attributed to the El-transition which we have not taken into account. 

For the correct description of the experimental data on the photodisintegration of the 
deuteron [40] (see also Chen and Savage [10]): 

ahD -> np) cxp (cu) = (1.300 ± 0.029) mb, 

a( 7 D^np) cxp M| w=2 76MeV = (1.474 ± 0.032) mb, 

a( 7 D^np) cxp H[ =445MeV = (2.430 ± 0.170) mb (5.5) 

one must include the contribution of the El-transition [10]. 

Since the cross section for the photodisintegration of the deuteron having been dis- 
cussed above does not contain the contribution of the El-transition, our result can be 
scarcely compared with the recent theoretical investigation of the photodisintegration of 
the deuteron carried out by Anisovich and Sadovnikova [41] within the dispersion relation 
approach based on the dispersion relation technique developed by Anisovich et al. [42]. 
This investigation represents a detailed analysis of the saturation of the cross section for 
the photodisintegration of the deuteron by different intermediate states valid mainly for 
the photon energy region uj > 50 MeV, where the contribution of the El-transition is 
important. 



6 The amplitude of the solar proton burning 

The analysis of the reaction p + p— >D + e + + z/ e within the NNJL model we should 

start with the evaluation of the effective Lagrangian C^~^ Be+l/a (x) of the transition p + p 

— > D + e + + u e . In the tree-meson approximation^] the effective Lagrangian C^Q~^ De " e (x) 

1 Below the analysis of weak nuclear reactions is carried out in the tree-meson approximation. The 
inclusion of chiral one-meson loop corrections goes beyond the scope of this paper and demands a separate 
publication. 
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is denned by the one-nucleon loop diagrams depicted in Fig. 3. The detailed evaluation 
of £g|^ De Uc (x) is given in Appendix. The result readsQ 

£ e P r Dc+ " (x) = -ig A C NN M N ^^ Dl(x) mxh 5 p(x)} [^(x)y*(l - T 5 )^*)]- (6.1) 

where G v = G F cos$ c with G F = 1.166 x HT 11 MeV~ 2 and $ c are the Fermi weak 
coupling constant and the Cabibbo angle cos^c = 0.975 

For the derivation of the effective Lagrangian Clg~^ De+Uc (x) we have used the effective 
Lagrangian £^f^ NN (x) responsible for low-energy transitions N + N — ► N + N defined 
by Eq. (li~T7T ). 

The matrix element of the transition p + p — > D + e + + u e we define by a usual way 

ln(U-\ +(U . ^,, ih \ I /'np-vDe+i/. 



d A x (D(kv)e + (k e+ )v e (k Ue )\£ir Uei ">(x)\p( Pl )p(p 2 )) = 

(o . , , . \ .M(n + p -»• D + e+ + u e ) 

= (2vr) S ( >{k D + k e + + k Vc -p x - p 2 ) , (6.2) 

^2E 1 V2E 2 V 2E B V 2E e+ V 2E Uc V 

where Ei(i = 1, 2, D, e + , u e ) are the energies of the protons, the deuteron, the positron 
and the neutrino, V is the normalization volume. 

The wave functions of the initial \p(pi)p(p 2 )) an d final (D{k- D )e + {k e +)v e {k Vc ) | state we 
take in the usual form 

|p(Pi)p(Pa)) = ^Op(Pi' t7 i) a p(^»' <7 ! l )l )' 
(D(k D )e + (k e +)v e (k Ve )\ = (0|a D (A; D , A D ) a e +(k e +, a) a Ve (k„ e , a'), (6.3) 

where a£(pi,Oi) and al{p 2 ,a 2 ) are the operators of creation of the protons. In turn, 

«d(^D) Ad), a e +(k e +, a e +) and a Ue (k Uc , a Uc ) are the operators of annihilation of the deuteron, 
the positron and the neutrino. 

The effective Lagrangian Eq. ( |6.1| ) defines the effective vertex of the low-energy nuclear 
transition p + p— >D + e + + z/ e 

UW(p + p^D + e + + z/ e ) = £ A M N C NN G v ^ <(A; D ) [u(k„ e )^(l --f 5 )v(k c+ )} 

x [u c (p 2 )-f 5 u(p 1 )}, (6.4) 

where e* (kv, Ad) is a 4-vector of a polarization of the deuteron, u(k Ue ), v(k e +), u{p 2 ) and 
u(pi) are the Dirac bispinors of the neutrino, the positron and the protons, respectively. 

For the evaluation of the matrix element Eq.( |6.4|) we have used the wave functions 
of the protons in the form of the plane waves. However, a real wave function of the pp 
pair in the 1 So state is defined by low-energy nuclear forces and Coloumb repulsion. In 
order to take into account both low-energy nuclear forces and Coulomb repulsion for the 
relative movement of the pp pair in the 1 So state we would sum up an infinite series of 
one-proton bubbles. The vertices of these one-nucleon bubbles are defined by 

V vp ^ w {k\ k) = C NN r pp (k' ) [u c ( P 2hMPi)} i> PP (k), (6.5) 



2 This result is obtained at zero contribution of the nucleon tensor current (see Appendix and discussion 
in the Conclusion). 
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where ip pp (k) and i/)* p (k') are the explicit Coulomb wave functions of the relative move- 
ment of the protons taken at zero relative distances, and k and k' are relative 3-momenta 
of the protons k = (p\ — p 2 )/2 and k' = (p[ — p 2 )/2 in the initial and final states. The 
explicit form of ?p PP (k) we take following Kong and Ravndal [33] (see also [43]) 

^Ak)=e-^ kr OT(l + -^-), (6.6) 



2kr c 

where 2rc = 2/M^a = 57.64 fm and a = 1/137 are the Bohr radius of the pp pair in the 
1 So state and the fine structure constant. The squared value of the modulo of ip pp (k) is 
given by 

|^ pW = C J W = ^_ 7 J_, ( 6 .7) 

where C (k) is the Gamow penetration factor [3,29,43]. 

By taking into account the contribution of the Coulomb wave function and summing 
up an infinite series of one-proton bubbles the expression Eq. ( |6.4| ) can be recast into the 
form 



[M c (p 2 hMpi)] % P (k) 



C N N f d 4 p Uj n ^,^ h|2 ^r„,5 1 5 1 



/ ^\^p(\p + Q\)\H^ 



(6.1 



16tt 2 J 7r 2 i I- M N -p-P-Q M N -p-Q- 



where P = Pi + p 2 = (2^Jk 2 + M^, ) is the 4-momentum of the pp-pair in the center 
of mass frame; Q = aP + bK = a (pi + p 2 ) + b (pi — p 2 ) is an arbitrary shift of virtual 
momentum with arbitrary parameters a and b, and in the center of mass frame K = 
Pi — P2 — (0,2k). The parameters a and b can be functions of k. 

The evaluation of the momentum integral we would carry out at leading order in 
the 1/Mn expansion or differently in the large Nq expansion [1] due to proportionality 
M N ~ N c valid in QCD with SU(N C ) gauge group at N c — > oo [22]. As a result we 
obtain 



^\^ p (\p + Q\)\ 2 k{j 5 



:7 



™ M N -p-P-Q' M N -j3-Q J 

= -8 a (a + 1) Ml + 8 (b 2 - a (a + 1)) k 2 - i 8vr M N k \^ pp (k) | 2 = 
= -8 a (a + 1) M 2 + 8 (b 2 - a (a + 1)) k 2 - i 8vrM N k C 2 (k). (6.9) 

Substituting Eq.( |6.9|) in Eq.( |6.8|) we get 

iM(p + p -> D + e + + u e ) = g A Mv C NN G Y 
x e;(k D ) [u(k Ua )r(l ~ 7>(M] m P2 h 5 u(Pi)] e-^/^c r(l + 

.GttNN „,2 , ^NN /,2 / , -i\\i2 • CnN^N 



2Ar c 



l-a(a + l)^M 2 +^f(^-a(a + l))^-^^f^A;C 2 (fc) . (6.10) 
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In order to reconcile the contribution of low-energy elastic pp scattering with low-energy 
nuclear phenomenology [43] we should make a few changes. For this aim we should rewrite 
Eq. (|6.10|) in more convenient form 



iM(p + p -> D + e H 



3<? 



V 



V 



4vr 2 



x e* (k D ) [u(k J>y»(l - 7 >(M] mp 2 )l h u{Pi)\ e 



ia (k) 



Co(fc) 



1 — a(a + 1) 



We have denoted 

e -n/Akr CT ^ 



2kr 



c 



(b 2 - a(a + l))k 2 -i 



M k ) C (k) , a (k) 



Cnn^n 
2tt 



arg 



r(i 



2kr 



c 



(6.11) 



(6.12) 



where a (k) is a pure Coulomb phase shift. 

Now, let us rewrite the denominator of the amplitude Eq. (|6.11 ) in the equivalent form 



| cos a (k) 



1 — a(a + 1 — - Mz, + 



NN 



sin cr (fc) N Cq(/c)| — 2 1 cos cr (A;) 



N 



2tt 



+ sin Co (A;) 1 — a(a + 1 



kC 2 (k) 



2tT 2 



2tt 2 



(6 2 



2tt 

a(a+l))k 2 ]} 



1 r 

2 



(6.13) 



where we have denoted 



1 

Z 



1 
2 



5 



C N nM n 2 
— sin do (re) — — — re C (re) 



+ cos Co (re) 1 - a(a + 1)— r-M^ + "■ 1,; 



2tt 2 



-i^*Cg(*) = cosa (re)^-^reC< 2 (re) 



2tt 2 

■»2 



2vr 

(fc 2 -a(a + l))re 5 



+ sin (To(fc) 



a a 



D— M N + — — (b -a(a + l))re 



2tt 2 



2vr 2 



(6.14) 



Here Z is a constant which we would remove by the renormalization of the wave functions 
of the protons, a e pp = (-7.8196 ± 0.0026) fm and r e pp = 2.790 ± 0.014 fm [44] are the S- 
wave scattering length and the effective range of pp scattering in the x So state with the 
Coulomb repulsion, and h(2krc) is defined by [43] 



h(2kr c ) = -7 + £n(2kr c ) + J2 



{n(l + 4n 2 re 2 r 2 7 ) ' 



(6.15) 



The validity of the relations Eq. ( |6.14| ) assumes the dependence of parameters a and b on 
the relative momentum k. 
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After the changes Eq.( 6.13p and Eq.( 6.14 ) the amplitude Eq.( |6.11 ) takes the form 



iM(p + p -> D + e + + z/ e ) = G v 9aM n C N n G v 



x 



e;(MK^7^(l-7>(M] 

(6.16) 



4vr 2 " 

Z [m c (p 2 )7 5 m(pi)]- 



r c 

Renormalizing the wave functions of the protons \J~Zu(j)<i) — > wfe) and \fZu(pi) — > 
we obtain the amplitude of the solar proton burning 

3flV 



iM(p + p -> D + e+ + i/ e ) = c/ A M N C NN G v 

C (A;) 



4vr 2 M 



X 



\ al p r; p k 2 + ^ h{2kr c ) + • a e pp fc C 2 (k) 



u'fofrVPi)]^ 2 ), (6.17) 



where F D (k 2 ) is given by Eq. fl2.16| ) and describes the spatial smearing of the deuteron 
coupled to the pp pair in the 1 So state. 

The real part of the denominator of the amplitude Eq. ( |6.17| ) is in complete agreement 
with a phenomenological relation [43] 

1 



ctgSlJk) 



C 2 (k) k 



pp 



(6.18) 



describing the phase shift S pp (k) of low-energy elastic pp scattering in terms of the S-wave 
scattering length a pp and the effective range r pp . Thus, we argue that the contribution of 
low-energy elastic pp scattering to the amplitude of the solar proton burning is described 
in agreement with low-energy nuclear phenomenology in terms of the S-wave scattering 
length a pp and the effective range r pp taken from the experimental data [44]. 



7 The astrophysical factor for the solar proton burn- 
ing 

The amplitude Eq. fl6.17|) squared, averaged over polarizations of the protons and 
summed over polarizations of final particles reads 

2 „2 li/rG ril 54Qe> z - 1 2/'j,2\ 



\M{v + P - D + e+ + v e )\* = G 2 y g{M^ C^ N F 2 (k 
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C 2 (k) 

X — 



1 - o <Mp* 2 + — h(2kr c ) + K p ) 2 k 2 C 4 (k) 



(E e+ E Ve -]-k e+ -k Ve ), (7.1) 



- n e r e k 2 4- ^ 

^Upp'ppft, T ■ W j , v - pp , 

where m e = 0.511 MeV is the positron mass, and we have used the relation g>y 
2rt 2 Q D M 2 . 

The cross section for the reaction p + p — > D + e + + u e is defined by 



(jPP^De+i 



^ = ^4^/ | - M(p + P ^ D + e+ + "« )|2 



x W < (fa + h - R - „) (2 ^ Eb ( ^ (2 g^ , (7.2) 
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where v is a relative velocity of the pp pair and = k c + + k Vc is a 4-momentum of the 
leptonic pair. 

The integration over the phase volume of the final De + z/ e state we perform in the 
non-relativistic limit 



d 3 ku d 3 k e 



il " U " (2n) 4 5^(k D + k e - Pl -p 2 ) (E c+ E Uc - \k c+ ■ k Uc ) 



r crfc D 

J (2vr) 3 2£ D '(2tt) 3 2E c+ (2tt) 3 2^ 

1 r w+T pp i (W + T ) 5 

where W = e D - (M n - M p ) = (2.225 - 1.293) MeV = 0.932 MeV and f = m e /(W + T pp ). 
The function /(£) is defined by the integral 



/(0 



30 



= 0.222 



(7.4) 



and normalized to unity at £ = 0. 

Thus, the cross section for the solar proton burning is given by 



+ -(T PP ) = f * a ^ff* 1 C;J N (W + T pp )= /( 



pp 



x 



PIV 



1 - « «; p C M n T pp + — H2r c ^M N T pp )f + {a c w fM^T vv C^^f x 



rc 



x 



1 = Spp(^pp) c -n/r C y/W\ 

l _ e -7r/r cv /M N T pp ^pp 



(7.5) 



The astrophysical factor S , PP (T PP ) reads 
bp P {T pp ) = a 



1280tt 4 



f( me ) 



x 



1-e "'"V'^i'i 1 
F 2 (M N T PP ) 



1 " £ «p P ^pM N T pp + ^ fc(2r c v/Mn^)] 2 + « p ) 2 M N T pp C 4 (^/M N T pp ) 



(7.6) 



At zero kinetic energy of the relative movement of the protons T pp = the astrophysical 
factor 5*^(0) is given by 



5pp(0) = ol 



9gjG^Q D M* 
1280tt 4 



C 2 N W 5 /(^) = 4.08 x l(T 25 MeVb. 



(7.7) 



The value Sp P (0) = 4.08 x 10 25 MeVb agrees well with the recommended value S'pp(O) = 
4.00 x 10~ 25 MeVb [29]. 
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Unlike the astrophysical factor obtained by Kamionkowski and Bahcall [30] the as- 
trophysical factor given by Eq. (\l.7\j does not depend explicitly on the S-wave scattering 
length of low-energy elastic pp scattering in the 1 So state. This is due to the normaliza- 
tion of the wave function of the pp pair. After the summation of an infinite series and by 
using the relation Eq. ( |6. 18|) we obtain the wave function of the pp pair in the l So state 



in the form 

% p (k)-e - a * p kC (ky (7 ' 8j 

that corresponds the normalization of the wave function of the relative movement of the 
pp pair used by Schiavilla et al. [31]. For the more detailed discussion of this problem we 
relegate readers to the paper by Schiavilla et al. [31] 0. 



8 The reaction z/ e + D^e +P + P 

The evaluation of the amplitude of the reaction v e + D — > e - + p + p is analogous 
to that of the amplitude of the solar proton burning. The result reads 

iM (u e + D - e" + p + p) = ffA M N | ^ C NN e*(Ato) [u(k e -)r{l - j 5 )u(k Ve )) 

x = [u(p 2 h 5 u c ( Pl )} F D (k 2 ). (8.1) 

1 " 5 + — h(2kr c ) +ia* p k C 2 (k) 

The amplitude Eq.fl3.ip squared, averaged over polarizations of the deuteron and summed 
over polarizations of the final particles reads 

\M(u e + D - e- + p + p)| 2 = S PP (Q) F ^ k 2 ) F + (Z, E e -) 



C 2 (k) 



< = e ~ uv ' v (E c -E v - -k e - • £„.)■ (8.2) 

1 " o«P fc2 + — mrc)] 2 + (al P fk 2 C*(k) 3 
Z r c 

where F + (Z,E e -) is the Fermi function [45] describing the Coulomb interaction of the 
electron with the nuclear system having a charge Z. In the case of the reaction u e + D 
— > e~ + p + p we have Z = 2. At a 2 Z 2 <C 1 the Fermi function F+(Z, E e ~) reads [45] 

+ F(Z 1 E c -)= 2 ^ e 2 , (8.3) 
1 _ e ^Ve- 



where i] e - = Za/v e - = ZaE c - / \J E 2 - — m 2 - and v e - is a velocity of the electron. 

For the evaluation of the r.h.s. of Eq.( |S.2|) we have also used the expression for the 
astrophysical factor 

Spp(0) = g^^ Cg N m«n^, (8.4) 



3 See the last paragraph of Sect. 3 and the first paragraph of Sect. 5 of Ref.[31]. 
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where tt De+Uc = (W/m e ff(m e /W) = 4.481 at W = 0.932 MeV. The function f(m e /W) 
is defined by Eq.flTl). 

In the rest frame of the deuteron the cross section for the reaction v e + D — > e~ + p 
+ p is defined by 



a 



i/ e D— >e pp 



\{2^5^\k T , + k Uc -p l 



X(z/ c + D^e~+p + p)| 2 

d 3 pi d 3 p 2 d 3 k e 



Vi ~ k. 



(2n) 3 2E 1 (2tt) 3 2£ 2 (2vr) 3 2£ c 



1.5) 



where E Uc , Ex, E 2 and E e - are the energies of the neutrino, the protons and the elec- 
tron. The integration over the phase volume of the (ppe~) state we perform in the 
non-relativistic limit and in the rest frame of the deuteron, 



d 3 pi d 3 p 2 



(27r) 4 6^\k D + k Vc - Pl -p 2 -k c 



2 J (2ir) 3 2E l (2tt) 3 2E 2 (2vr) 3 2£ e - 

C 2 (Jm n T pp ) F*(M n T pp ) F + (Z, E 



rc 



1 " oWn^pp + h(2r c JM N T pp ) + (a*) 2 M N T PP C^(J M N T PP ) 



dT e -dT pp 5(E Vc - E th - T c - - T pp )^fT e -T pp l 1 + -2- 1 Jl + ^ 



E 2 
T 

1 n 



C*UM N T PP ) i^(M N T pp ) F+(Z, E e -] 



rc 

X 7 /2 /«, \3/2 



1 " 2 a PP r PP M N T pp + 7^ h(2r c ^/M N T pp ) + (o^MnTppCoVMnTpp) 
£, c M 3 (E lh \ U2n 



128tt 3 M N 



tb 



(y-1) fippe-(y), 



5.6) 



where T e - is the kinetic energy of the electron, E t h is the neutrino energy threshold: 
E th = e D + m e - (M n - M p ) = (2.225 + 0.511 - 1.293) MeV = 1.443 MeV. The function 
Q ppe -(y), where y = E Uc /E th and E th = e D — (M n — M p ) +m e = 1.443 MeV, is determined 
by the integral 



fi P pc- (y) = / dxy/x(l-x)(l + ^(y - 1)(1 - *)) J 1 + ^L(y - 1)(1 - 



C^M N E th (y - 1) rc) F£(M N E th (y - 1) s) m c + £ th (y - 1) (1 - x)) 

1 " .<p<p M N^th (y - 1) x + h(2r c JM N E th (y-l)x 



rc 



+ (a e pp ) 2 M N £ th (y - 1) xC*(JM N E th (y - 1) x) 



1.7) 



where we have changed the variable T pp = (E Ve — E t ^) x. 
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The cross section for the reaction z/ e + D — > e +p + pis defined by 








= 4 MeV) 




—pp^ 


= 6 MeV) 






= 8 MeV) 




C_PP (^ C = 


= 10 MeV) 



= 3.69 x 10 5 S pp (0) (y - l) 2 fi ppe - (y), (8.8) 

where 5 PP (0) is measured in MeV cm 2 . For S pp (0) = 4.08 x 10" 49 MeVcm 2 Eq.(|7\7|) the 
cross section ( j I ' eI - ) ~ >e pp^g^j reads 

a ^- P P {E j = 1.50(y-l) 2 n ppc -(y) 10- 43 cm 2 . (8.9) 

Recently the calculation of the cross section for the reaction z/ e + D^e _ +p + p has 
been carried in Ref. [46] within the PMA and tabulated for the neutrino energies ranging 
over the region from threshold up to 170 MeV. Since our result is valid for much lower 
neutrino energies, we give the numerical values of the cross section only for the energies 
from the region 4 MeV < E Uc < 10 MeV: 

= 2.46(1.577) x 10~ 43 cm 2 , 
= 10.04(6.239) x 10" 43 cm 2 , 
= 2.37(1.463) x 10" 42 cm 2 , 
a ».u-«rw(E Va = 10 MeV) = 4.39(2.708) x 10" 43 cm 2 . (8.10) 

The data in parentheses are taken from Table 1 of Ref. [46]. Thus, on the average the 
cross section a u ^^ c pp (E Uc ) calculated in the NNJL model by a factor of 1.6 is larger 
compared with the PMA ones. 

Since the amplitude of the reaction u e + D — * e~ + p + p is completely defined by 
the amplitude of the solar proton burning p + p^D + e + + z/ c that is described in 
the NNJL model in agreement with other theoretical approaches, our prediction for the 
cross section for the neutrino disintegration of the deuteron Eq. (|8.10|) is a challenge to 
the experiments planned at SNO [35]. 

In order to compare the cross section for the neutrino disintegration of the deuteron u e 
+ D^e~+p + p with experimental data planned to be obtained by SNO Collaboration 
we should average it over the 8 B solar neutrino energy spectrum produced by the (3 decay 
8 B — > 8 Be* + e + + u e in the solar core. Using the 8 B solar neutrino energy spectrum [47] 
and integrating over the region 5 MeV < E Uc < 15 MeV [48] , where the lower bound is 
related to the experimental threshold of experiments at SNO and the upper one is defined 
by the kinematics of the decay 8 B — ► 8 Be* + e + + u e , we get 

(a u ^ c ~ pp (E Uc )) H s B) = 2.62 x 10" 42 cm 2 . (8.11) 

For the comparison of the theoretical cross section with the experimental one measured 
at SNO one should take into account a possible decrease of the experimental value in the 
case of the existence of neutrino flavour oscillations [27,48]. 



9 The astrophysical factor for the pep process 

In the NNJL model the amplitude of the reaction p + e~ + p — > D + u e or the 
pep-process is related to the amplitude of the solar proton burning p + p^D + e + + 
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v e Eg . ( |6. 17] ) as well as to the amplitude of the neutrino disintegration of the deuteron u e 
+ D — > e~ + p + p. Indeed, the effective Lagrangians 

^eff D ^ e PP (%) are defined by the anomaly of the one-nucleon loop triangle AAV-di&grams 
with two axial-vector (A) and one vector (V) vertices (see Appendix). The amplitude of 
the pep-process reads 

iM(p + e" + p - D + u e ) = g A M N C nn G Y ^ e ;(k B ) [u(k Ue )^{l - ^)u(k e -)} 



x Mp 2 hMPi)] Mk 2 )- (9.1) 

1 - 5 « p £; 2 + -2£ fc(2*rc) + • a c pp k C 2 (k) 

Z Tq 



The amplitude Eq.( |9.1| ) squared, averaged and summed over polarizations of the inter- 
acting particles is defined by 



\M( V + e- + p - D + v e )\> = glM* C 2 N G% ^ F 2 (k 2 ) F + (Z, E e -) 
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x = Q c ^ (E c+ E Uc - \k e+ ■ Q, (9.2) 

1 " 5 alAk 2 + h(2kr c )} + (al P ) 2 k 2 C^k) 

Z Tq 

where F + (Z, E e -) is the Fermi function given by Eq.([S.3|). 

At low energies the cross section cr pc p ^ Duc {T pp ) for the pep-process can be determined 
as follows [49] 

(T PP ) = ij^j / s n ( k r ) J | M( P T e- + p = D + V , W 

{2mB + ^ _ „ _ p2 _ (9 . 3 ) 

where (7 = 2 is the number of the electron spin states and v is a relative velocity of the 
pp pair. The electron distribution function n(k e -) can be taken in the form [45] 

n(jfe e -) = e v-T e -/kT c ^ ( 9 4 ) 

where fc = 8.617 x 10~ n MeV K -1 , T c is a temperature of the core of the Sun. The 
distribution function n(k e ~) is normalized by the condition 

9 jy~j3 n &-) = n e-, (9.5) 

where n e - is the electron number density. From the normalization condition Eq.( |9.5|) we 
derive 

p 47r 3 n e - 
~ (2tt m e kT c f' 2 ' 

The astrophysical factor S , pep (0) is then defined by 



S pep (0) = S pp (0) £ ^— ^ / c*V e- T *-/ kT cF + (Z,E e -). (9.7) 



27T fi D e+i/ e 
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For the ratio S pep (0) / S pp (0) we obtain 

ff pep (0) _ 2 3 /y/ 2 ( aZn e - \ (E th V [^ r ( 7 ^A (qR) 
S PP (0) /pp(0) ^ ml J \m c ) \jkT c \ V kTj' 1 ' } 

We have set / pp (0) = fiDe+j/ c /30 = 0.149 [45] and the function I(x) introduced by Bahcall 
and May [45] reads 




The relation between the astrophysical factors ^^(0) and 5 PP (0) given by Eq. (|9~8| ) is in 
complete agreement with that obtained by Bahcall and May [45]. 

By virtue of the direct relation between the amplitudes of the pep-process and the 
reaction for the neutrino disintegration of the deuteron u c + D — >e~+p + p that we 
have in the NNJL model the agreement with the result obtained by Bahcall and May [45] 
is on favour of our predictions for the cross section for the reaction u e + D — > e~ + p + 

P- 



10 The reaction z/ e + D^e + + n + n 

The effective Lagrangian £g| D ^ e+nn (x) of the low-energy nuclear transition z/ e + D 
— > e + + n + n we evaluate by analogy with C^~* Uc (x) (see Eq.(A.27) of Appendix) 
through the one-nucleon loop exchanges and at leading order in the large Nc expansion. 
The effective Lagrangian £g| D ^ e nn (x) is defined by the anomaly of the one-nucleon loop 
triangle AA^Z-diagram as well as the effective Lagrangian C^^ Dc Uc (x). The result reads 

£^ c+nn (x) = igA M^C NN ^^ D,(x) [n(x) 7 V(x)] [^.(x)y(l - 7 5 )<M*)]. (10.1) 

The amplitude of the reaction z/ e + D — > e + + n + n we obtain in the form 

G 3 1 

tM(u c + B + n + n) = _^ A M N C NN ^^ , F D {k 2 ) 

V2 2vr^ 1 _ £ p, jn u 

^ 2 nn u, nn fu i u "nil 

xe,(k D ) [v(k Pc )Y(l - l>(k c+ )} [u(p 2 ) 7 V( Pl )], (10.2) 

where the form factor Fj)(k 2 ) describes a spatial smearing of the deuteron Eq. (|2.16| ). The 
factor 1/(1 — \r nn a nn k 2 + i a nn k) gives the contribution of low-energy nuclear forces to the 
relative movement of the nn pair in the 1 So state. The result is obtained by the summation 
of an infinite series of the one-neutron bubbles evaluated at leading order in the large 
Nc expansion. Since we work in the isotopical limit, we set a nn = a np = — 23.75 fm and 
r nn = r np = 2.75 fm. The recent experimental values of the S-wave scattering length and 
the effective range of low-energy elastic nn scattering are equal to a nn = (— 18.8 ± 0.3) fm 
and r nn = (2.75 ±0.11) fm [50,51]. 
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The amplitude Eq. (|10.2|) , squared, averaged over polarizations of the deuteron and 
summed over polarizations of the final particles, reads 

\M[y e + D^e+ + n + n)| 2 = — j [E e +E Dc - -k c + ■ k Pc ). (10.3) 

(l 2^ nnCtnn ^ ) a nn^ 

The cross section for the reaction u c + D — > e + + n + n is defined by 

1 



^ C nn (£pJ = ^7 ||^c + D-eHn + n) 

-(2n) 6"{Q + k Pe - Pl -p 2 - k e+ ) (27r)32Ei (27r)32E2 (27r)32Ee+ , (10.4) 

where E D , E Pe , Ei, E 2 and E e + are the energies of the deuteron, the anti-neutrino, the 
neutrons and the positron. The integration over the phase volume of the (nne + ) state we 
perform in the non-relativistic limit and in the rest frame of the deuteron 

1 r d 3 p x d 3 p 2 d 3 k c + \4 X u) fn . u / \ 

2 J (2n) 3 2 El (2n) 3 2E 2 (2n) 3 2E c+ ^ 5 {Q + ^ ~ Pl ~ P2 " h+) 

~ \k + ■ Q F 2 (M N T nn ) Ep M3 / Eth y/ 2 /2m c \ 3/2 o 
'l l r n MT\ 2 + n 2 MT ~ 1024tt 2 I^mJ \E th ) 



dT r +dT, 



T e +T nn Fp (M N T n 

'"' '1 \ 2 2 

1 - -^nna n nM N T nn + GrM N T nn 



// 



( 1+ i:)\^£^-^- T = + - r » 



E V M 3 {E th \ 7/2 {2m e \ 3/2 



1024tt 2 Im n J i 



th 



(j/-!) 2 ^^), (10-5) 



where T nn is the kinetic energy of the nn pair in the 1 So state, T e + and m e = 0.511 MeV are 
the kinetic energy and the mass of the positron, y = E Pc /E t h and E t h is the anti-neutrino 
energy threshold of the reaction u e + D — > e + + n + n: E th = e D + m e + (M n — M p ) = 
(2.225 + 0.511 + 1.293) MeV = 4.029 MeV. The function tt QnPc (y) is defined by 



8 r 

tt an p B (y) = - J dx y 



x (1 - x) F 2 (M N £ th (y-l)x) 



1 - ^r nn a nn M N E th (y -l)xj + a 2 n M N E th (y - 1) x 

x (l + ^(y - 1)(1 - x)) ^l + ^-lXl-x), (10.6) 

where we have changed the variable T nn = (£ , Pc — Z? t h) x. The function /(i/) is normalized 
to unity at y = 1, i.e. at threshold -E Pc = £? t h- Thus, the cross section for the reaction u e 
+ D— >e + + n + n reads 



^D^e+nn^j = ^ (y _ ^ (10 . 7) 



29 



where a amounts to 

a -g D C NN 5127r4 y— J ^—j -4.58x10 cm. (10.8) 

The experimental data on the anti-neutrino disintegration of the deuteron are given in 
terms of the cross section averaged over the anti-neutrino energy spectrum [36] : 

(a PcD ^ c+nn (E Pc )) cxp = (9.83 ± 2.04) x 10" 45 cm 2 . (10.9) 

In order to average the theoretical cross section Eq. ( |10.7| ) over the anti-neutrino spectrum 
we should use the spectrum given by Table YII of Ref . [36] . This yields 

( a *cD— + ™(£_J) = n x 10 -45 cm 2_ ( 1(U0 ) 

The theoretical value Eq.( |10.10"D 4 agrees well with the experimental one Eq.( |10.9[ ). 



11 The reactions v e (v e ) + D —> v e (v e ) + n + P 

The amplitude of the neutrino disintegration of the deuteron caused by neutral weak 
current i/ e + D^z/ e + n + p evaluated through one-nucleon loop exchanges and at 
leading order in the large Nc expansion reads (see Eq.(A.30) of Appendix): 

G 3 1 

iM{v e + D -> u e + n + p) = - g A M N C NN -| ^ j F D {k 2 ) 

1 2^np'-^np^' i ^np k 

x e,(k D ) [u(k')Y(l ~ ^)u(K e )MP2h 5 u c (pi)\. (11.1) 



The amplitude Eq.( ll.lf ) squared, averaged over polarizations of the deuteron, summed 



over polarizations of the nucleons reads 

|A4(z/ e + D -»• z/ e + n + p)| 2 = — j — -—[E v E Ve - -k Ve ■ k Uc j . (11.2) 

^1 2^ n P an P^ J a np^ 

In the rest frame of the deuteron the cross section for the reaction z/ e + D^z/ c + n + p 
is defined by 



y e D^i/ c iip 



(2vr) 4 5( 4 )(A; D + k Vc - Pl -p 2 -k' v ) , ^f 1 . ^f 2 . ^ . (11.3) 
v 7 v c ^ ^ ^ (27r) 3 2Ei (2tt) 3 2£ 2 (2tt) 3 2^ c v ; 

The integration over the phase volume of the (npz/ e ) state we perform in the non- 
relativistic limit and in the rest frame of the deuteron, 

r (P Pl d 3 p 2 d 3 k' Uc 4 (4) 



4 This result is obtained at zero contribution of the nucleon tensor current (see Appendix and the 
discussion in the Conclusion). 
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2 ' np u, np-"'- 1 lM J- npy i "up . . 

f (tf fa " 1)7/2 n - iv) - (1L4) 

The function f2 npi , c (?/), where y = E Pc /E th and E th = e D = 2.225 MeV is threshold of the 
reaction, is defined by the integral 

( , 105 r y/j(l-x) 2 1 m 

tonpuAy) = -T7T / dx ; : — ~% jvtt, — TT^i> ( 1L5 ) 



!6 J (. lr np a np ,2 a (l + (y-l)x) 2 ' 

'1-- — 2— (y-l)x 



Tf 



D ' D 

where we have changed the variable T np = (E Pc — E t h) x and used the relation M-^E t h = 
1/^d at Eth = £d- The function Q npi/e (y) is normalized to unity at y = 1, i.e. at threshold 
Ep c = E t h- 

The cross section for the reaction z/ e + D— >z/ e + n + p reads 

a u ^ v ^{E Ve ) =a (y- I) 7 / 2 fl^y), (H-6) 

where <r amounts to 

, = g DC J N M^lgJ /2 = 1 . 8 4xlO-c^ (11.7) 

In our approach the cross section for the disintegration of the deuteron by neutrinos u e + 
D — ► z/ e + n + p coincides with the cross section for the disintegration of the deuteron by 
anti-neutrinos u e + D — > u e + n + p, a UcD ^ Uclip (E Uc ) = a PcD ^ Pcnp (E Pc ). Therefore, we have 
an opportunity to compare our result with the experimental data on the disintegration 
of the deuteron by anti-neutrinos [36]. The experimental value of the cross section for 
the anti-neutrino disintegration of the deuteron z/ e + D — >z/ e + n + p averaged over the 
anti-neutrino spectrum reads [36]: 

(a PcD ^ Pcnp (£ p J) cxp = (6.08 ± 0.77) x 10" 45 cm 2 . (11.8) 

By using the anti-neutrino spectrum given by Table YII of Ref. [36] for the calculation of 
the average value of the theoretical cross section Eq.( 11.6 ) we obtain 

(a PcD ^ np {E Pc )) = 6.28 x 10~ 45 cm 2 . (11.9) 

The theoretical value Eq.( |11.9[) p| agrees well with the experimental one Eq. (|11.8|) . 



The cross section for the neutrino disintegration of the deuteron u e + D — > v c + n 
+ p averaged over the 8 B solar neutrino spectrum [47] for energy region 5 MeV < E Ve < 
15 MeV is given by 

(a UeB ^ np (E Ve )) H s B) = 1.85 x 10~ 43 cm 2 . (11.10) 

This result can be directly compared with the experimental data obtained by SNO, since 
the averaged value for the cross section for the reaction z/ e + D— >z/ e + n + p caused by 
the neutral weak current should not depend on whether neutrino flavours oscillate or not 
[27,48]. Of course, if none sterile neutrinos exist [27,48]. 

5 This result is obtained at zero contribution of the nucleon tensor current (see Appendix and the 
discussion in the Conclusion). 
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12 Conclusion 



We have considered the description of a dynamics of low-energy nuclear forces within 
the Nambu-Jona-Lasinio model of light nuclei (the NNJL model) for low-energy electro- 
magnetic and weak nuclear reactions with the deuteron. We have shown that the NNJL 
model enables to describe in a reasonable agreement with both experimental data and 
other theoretical approaches all variety of low-energy electromagnetic and weak nuclear 
reactions with the deuteron in the final or initial state coupled to a nucleon-nucleon (NN) 
pair in the state. In the bulk the reaction rates for the neutron-proton radiative 
capture n + p — > D + 7 for thermal neutrons, the photomagnetic disintegration of the 
deuteron 7 + D — > n + p, the solar proton burning p + p — > D + e + + z/ e , the pep-process 
p + e~ + p — > D + u e , the disintegration of the deuteron by neutrinos and anti-neutrinos 
caused by charged z/ e + D — >e - +p + p and z/ c + D^e + + n + n and neutral v e (v e ) 
+ D — > v e (v e ) + n + p weak currents are calculated in agreement with other theorical 
approaches and experimental data. 

When matching our results with those obtained in the PMA we would like to emphasize 
a much more simple description of the NN interaction responsible for the transition N + 
N — > N + N and a substantial simplification of the evaluation of matrix elements of low- 
energy nuclear transitions near thresholds of the reactions where the NN pair in the 1 So 
state is produced or absorbed with a relative momentum p comparable with zero. Such a 
simplification is rather clear in the NNJL model where with a good accuracy the deuteron 
can be considered within a quantum field theoretic approach as a point-like particle. 
Indeed, the spatial region of the localization of the NN pair is of order of 0(1 /p). Near 
thresholds the effective radius of the deuteron r D = 4.319 MeV is much smaller than 1/p, 
t-q <C 1/p. This yields that the NN pair does not feel the spatial smearing of the deuteron 
and couples to the deuteron as to a point-like particle. A correct description of strong 
interactions of the point-like deuteron coupled to the NN pair in the 1 So state is guaranteed 
then by the one-nucleon loop exchanges with dominant contributions of the nucleon- 
loop anomalies. This implies that the effective Lagrangians C^ T>1 {x), Cl^ Be+Uc (x), 
£ef^ C ~ PP (%), £eff~ P " D " e (z), C p e f^ e+nn (x) and C" e f-* Penp (x) are well defined. Thus, the 
procedure of the derivation of effective Lagrangians of low-energy nuclear transitions in 
the NNJL model treating the deuteron as a point-like particle coupled to nucleons and 
other particles through one-nucleon loop exchanges seems to be good established. We 
argue that the application of this procedure should get correct results for the derivation 
of effective Lagrangians of any other low-energy nuclear transitions, effective vertices, of 
nuclear reactions. 

Some problems occur for the evaluation of the amplitudes of nuclear reactions de- 
manding the continuation of matrix elements of low-energy nuclear transitions defined by 
the effective Lagrangians to the energy region far from thresholds. 

The continuation of matrix elements of low-energy nuclear transitions demands in 
the NNJL model: 1) the spatial smearing of the NN pair in the 1 So state and 2) the 
spatial smearing of the deuteron caused by the finite value of the effective radius r^. The 
spatial smearing of the NN pair in the 1 So state can be carried out by a summation of 
an infinite series of one-nucleon bubbles describing rescattering or differently a relative 
movement of the NN pair either in an initial or a final state of a nuclear reaction. The 
result of the NN rescattering can be expressed in terms of S-wave scattering lengths and 
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effective ranges in complete agreement with nuclear phenomenology of low-energy elastic 
NN scattering in the 1 So state. However, for the description of the spatial smearing of 
the deuteron the abilities of the NNJL model are rescticted and most what one can do at 
present level of the development of the model is to introduce the spatial smearing of the 
deuteron phenomenologically in the form of the Fourier transform of the approximate 3 Si 
wave function of the deuteron normalized to unity at p = 0: F-£>(p 2 ) = 1/(1 + r r>P 2 )- We 
have chosen a simplest form of the spatial smearing of the deuteron. Of course, F^(p 2 ) 
can be taken in the more complicated form of the Fourier transform of the explicit wave 
function of the deuteron. Of course, such a dependence is not absolute and the spatial 
smearing of the deuteron can be taken into account in the form of phenomenological form 
factors as it has been done by Mintz [52], for example. 

We would like to emphasize that the cross section for the Ml-capture n + p — > D + 
7 for thermal neutrons has been calculated by accounting for the contributions of chiral 
one-meson loop corrections and the A(1232) resonance. The total cross section for the 
Ml-capture has been found dependent on the parameter Z defining the 7rNA coupling 
constant off-mass shell of the A(1232) resonance: er(np — > D7)(T n ) = 325.5 (1 + 0.246 (1 — 
2Z)) 2 mb Eq. (|4.15|) . In order to fit the experimental value of the cross section we should 
set Z = 0.473. This agrees with the experimental bound \Z\ < 1/2 [18]. At Z = 1/2 
the contribution of the A(1232) resonance to the amplitude of the Ml-capture is defined 
only by the nucleon tensor current. Setting Z — 1/2 that is favoured theoretically [15] 
we calculate the cross section for the Ml-capture o~(np — > D7)(T n ) = 325.5 mb agreeing 
with the experimental data with accuracy better than 3%. 

When matching our result for the cross section for the Ml-capture with the recent 
one obtained in the EFT approach by Chen, Rupak and Savage [9]: cr(np — > Dj)(T n ) = 
(287.1 + 6.51 9 Xi) mb (see Eq.(3.49) of Ref. [9]), we accentuate the dependence of the cross 
section on the parameter *L\ undefined in the approach. This parameter has to be fixed 
from the experimental data [9]. In the NNJL model the parameter *L\ can be expressed 
in terms of the Z parameter as follows: = 5.90 + 24.60 (l-2Z)+3.03 (1-2Z) 2 . Thus, 
in the NNJL model the parameter acquires a distinct meaning of the contribution of 
the A(1232) resonance. 

The obtained result for the Ml-capture n + p^D+7we have applied to the 
analysis of the photomagnetic disintegration of the deuteron 7 + D — ► n + p. Due 
to time-reversal invariance the cross section for the photomagnetic disintegration of the 
deuteron can be directly expressed through the cross section for the Ml-capture. We 
have compared the numerical values of the cross section a^D — ► np)(u;) calculated in 
the NNJL model with the results obtained by Chen and Savage [10] and found a good 
agreement. Nevertheless, it should be emphasized that in the critical region of the photon 
energies u < 2 Ed = 4.45 MeV restricting the energy region of the dominance of the 
photomagnetic disintegration of the deuteron the cross section calculated in the NNJL 
model falls steeper than the cross section obtained in the EFT. However, in this region 
the dominant role should be attributed to the El-transition [10], which we have not 
considered. The comparison of our results on the photomagnetic disintegration of the 
deuteron with the experimental data demands the inclusion of the El-transition as well. 
This is planned in our forthcoming publications. 

The effective coupling constants of low-energy weak transitions p + p^D + e + + 
u e , p + e~ + p — > D + u e , u e + D — > e~ + p + p, u e + D — ► e + + n + n and u e (u e ) + 
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D — > v c {v c ) + n + p have been found dependent on an arbitrary parameter £ in the form 
9v — ¥ 9v (1 + caused by the contribution of the nucleon tensor current [1]. 

Since at £ = we get the value of the astrophysical factor S pp (0) = 4.08 x 10 -25 MeVb 
in complete agreement with the recommended one S pp (0) = 4.00 x 10~ 25 MeVb related to 
the Standard Solar Model [27,28], any non-zero value of £ should lead to an Alternative 
Solar Model (ASM). The problem of the formulation of the ASM goes beyond the scope 
of this paper. Therefore, in order to dwell within the Standard Solar Model [27,28] we 
have to set £ = 0. 

Setting £ = we have shown that all low-energy weak nuclear reactions of astrophysi- 
cal interest (i) the solar proton burning p + p^D + e + + z/ e , (ii) the pep-process p + e~ 
+ p — > D + v c and (iii) the disintegration of the deuteron by neutrinos and anti-neutrinos 
caused by charged u c + D — >e~+p + p and u c + D — >e + + n + n and neutral u c {u c ) + 
D — > u c {u c ) + n + p weak currents are described in the bulk in agreement with other the- 
orical approaches and experimental data. The effective Lagrangians of low-energy weak 
nuclear transitions p + p — > D + e + + u c , p + e~ + p — > D + u e , v e + D — > e~ + p + p, 
z/ c + D— >e + + n + n and u e (i/ e ) + D — > v c {v c ) + n + p are determined by the anomaly 
of the one-nucleon loop triangle AAV-diagrams. This confirms the statement argued in 
the NNJL model concerning the dominant role of the one-nucleon loop anomalies for the 
description of low-energy nuclear forces within a quantum field theoretic approach. 

We have shown that the contributions of low-energy elastic pp scattering in the 1 So 
state with the Coulomb repulsion to the amplitudes of the reactions p + p^D + e + + 
v e , v G + D — > e~ + p + p and p + e~+p^D + z/ c are described in the NNJL model in 
full agreement with low-energy nuclear phenomenology in terms of the S-wave scattering 
length and the effective range. The amplitude of low-energy elastic pp scattering has been 
obtained by summing up an infinite series of one-proton loop diagrams and evaluating the 
result of the summation at leading order in the large Nq expansion. The same method 
has been applied to the evaluation of the contribution of low-energy nuclear forces to the 
relative movements of the nn and np pairs, respectively, for the reactions v c + D — > e + 
+ n + n and v c (u c ) + D — > v c (u c ) + n + p. This has given the amplitudes of low-energy 
elastic nn and np scattering described in terms of S-wave scattering lengths and effective 
ranges in agreement with low-energy nuclear phenomenology [26] as well. 

The astrophysical factor S , pep (0) for the pep-process, p + e~ + p — > D + u c , evaluated 
relative to 5 , pp (0) is found in full agreement with the result obtained by Bahcall and May 
[45]. 

The cross sections for the anti-neutrino disintegration of the deuteron caused by 
charged u e + D^e + + n + n and neutral P e + D— >P e + n + p weak currents 
and averaged over the anti-neutrino spectrum (a UcD ~* c+nn (E Pc )} = 11.56 x 10~ 45 cm 2 and 
(a UcB ^ Ucnp (Ep c )) = 6.28 x 10~ 45 cm 2 agree well with recent experimental data 

{a 9 ^ e+nn (E Pc )} cxp = (9.83 ± 2.04) x 10~ 45 cm 2 , 
(<7 PeD -^ np (£pJ) exp = (6.08 ± 0.77) x 10~ 45 cm 2 

obtained by the Reines's experimental group [36]. 

The cross sections for the reactions z/ c + D-^e + + n + n and z/ e + D^z/ e + n + p 
have been recently calculated by Butler and Chen [53] in the EFT approach. The obtained 
results have been written in the following general form a — (a + bL 1; \) x 10~ 42 cm 2 (see 
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Table I of Ref. [53]), where a and b are the parameters which have been calculated in 
the approach, whereas L 1; a is a free one. In the NNJL model the appearance of the free 
parameter is related to the contribution of the nucleon tensor current [1] that effectively 
leads to the change of the coupling constant gy — > g\ (1 + £), where £ is an arbitrary 
parameter. The best agreement with the recommended value of the astrophysical factor 
for the solar proton burning [29] and the contemporary experimental data [36] on the 
cross sections for the anti-neutrino disintegration of the deuteron z/ c + D^e + + n + n 
and z/ e + D^z/ c + n + p caused by charged and neutral weak current, respectively, we 
obtain at £ = (see Appendix). 

We would like to emphasize that the contribution of the A(1232) resonance to the 
amplitudes of the low-energy weak nuclear reactions p + p^D + e + + z/ c ,p + e~ 
+ p^z/ c + D, z/ c + D-^e~+p + p, z/ c + D— >e + +n + n and u c + D — > u c 
+ n + p can be neglected. In fact, the contribution of the A(1232) resonance to the 
amplitudes of these reactions is of order of the momentum of the leptonic pair. This is 
due to the gauge invariance of the effective interactions ANW + and ANZ which should 
be proportional to = d^Wj — d u W* and = d il Z v — d u Z^, respectively. Since the 
amplitudes of the low-energy weak nuclear reactions under consideration are defined by 
the Gamow-Teller transitions, the terms proportional to the momentum of the leptonic 
pair give negligible small contributions. Thus, the cross sections for low-energy weak 
nuclear reactions enumerated above do not depend practically on the uncertainties of the 
parameter Z. 

This distinguishes low-energy weak nuclear reactions with the deuteron from the 
neutron-proton radiative capture analysed in the NNJL model. Unlike the low-energy 
weak nuclear reactions the contribution of the A(1232) resonance to the amplitude of 
the neutron-proton radiative capture is essential for the explanation of the experimental 
data. 

The cross section for the reaction z/ e + D^e~+p + p has been evaluated with 
respect to S , pp (0). We have found an enhancement of the cross section by a factor of 1.6 
on the average relative to the results obtained in the PMA. It would be important to 
verify this result for the reaction v c + D — > e - + p + p in solar neutrino experiments 
planned by SNO. Indeed, first, this should provide an experimental study of S'pp(O) and, 
second, the cross sections for the anti-neutrino disintegration of the deuteron caused by 
charged z/ c + D^e + + n + n and neutral z/ e + D^z/ c + n + p weak currents have 
been found in good agreement with recent experimental data [36]. 

For the comparison of the cross sections for the reactions z/ e + D — > e _ + p + p and 
z/ c + D^z/ c + n + p caused by charged and neutral weak currents, respectively, we have 
averaged the cross sections over the 8 B solar neutrino energy spectrum at energy region 
5 MeV < E Ve < 15 MeV, where the lower bound is related to the experimental threshold 
of the experiments at SNO and the upper one is defined by the kinematics of the (3 decay 
8 B — > 8 Be* + e + + v e being the source of the 8 B neutrinos in the solar core. We have 
obtained 

( a -^-pp(E Uc )) H s B) = 2.62 x l(T 42 cm 2 , 
K cD ^ cnp (^J) <E ,(B B ) = 1-85 x l(T 43 cm 2 . 

The experimental value for the cross section for the reaction z/ e + D — >e~+p + p caused 
by the charged weak current can, in principle, differ from the theoretical one due to a 
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possible contribution of the neutrino flavour oscillations [27,48]. In turn, the averaged 
value of the cross section for the reaction z/ c + D^z/ c + n + p caused by the neutral 
weak current should be directly compared with the experimental data, since it should not 
depend on whether neutrino flavours oscillate or not. Of course, the former is valid only if 
there is no so-called sterile neutrino [27,48] having no interactions with Standard Model 
particles [39]. 

Concluding the paper we would like to emphasize that the NNJL model conveys the 
idea of a dominant role of one-fermion loop (one-nucleon loop) anomalies from elementary 
particle physics to the nuclear one. This is a new approach to the description of low- 
energy nuclear forces in physics of light nuclei. In spite of almost 30 years have passed 
after the discovery of one-fermion loop anomalies and application of these anomalies to 
the evaluation of effective chiral Lagrangians of low-energy interactions of hadrons, in 
nuclear physics fermion-loop anomalies have not been applied to the analysis of low- 
energy nuclear interactions and properties of light nuclei. The contributions of one- 
nucleon loop anomalies are strongly related to high-energy NN fluctuations of virtual 
nucleon fields [1,54]. An important role of NN fluctuations for the correct description 
of low-energy properties of finite nuclei has been understood in Ref . [55] . Moreover, NN 
fluctuations have been described in terms of one-nucleon loop diagrams within quantum 
field theoretic approaches, but the contributions of one-nucleon loop anomalies have not 
been considered. The NNJL model allows to fill this blank. Within the framework of the 
NNJL model we aim to understand, in principle, the role of nucleon-loop anomalies for 
the description of a dynamics of low-energy nuclear forces at the quantum field theoretic 
level. 
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Appendix. The effective Lagrangian of the transition 

p + p^D + e+ + z/ e 

The reaction p + p^D + e + + z/ e runs through the intermediate W-boson exchange, 
p + p — > D + W + — > D + e + + u e . In the NNJL model we determine this transition in 
terms of the following effective interactions [1,39] 



•-npD 



[x) = —igy\p c (x)^n(x) — n c (x)^p(x)} D* (x), 



:iy pp (x) = \ C NN {[p(x) ^Y(x)] mx) ^ 5 p(x)\ + \p(x) 7 V (x)} \p^x) ^p(x)]}, 

9w 



C npW (x) = costfc [n(x)Y(l ~ 9Al 5 )p(x)} W~(x). (A.l) 

The transition W + — * e + + u e is defined by the Lagrangian [39] 

£, c e+w(*) = -ff= [^ e (x)r(l ~ W+(x). (A.2) 

The electroweak coupling constant gw is connected with the Fermi weak constant Gp and 
the mass of the W boson M w through the relation 



<?w _G F (A3) 



8M^ y/2' 

For the evaluation of the effective Lagrangian C P g~^ Dc+Uc (x) it is convenient to use the 
interaction 

£n P w(z) = [n(x)Y(l - 9Al 5 )p(x)} W~(x), (A.4) 

and for the description of the subsequent weak transition W + — > e + + v c to replace the 
operator of the W-boson field by the operator of the leptonic weak current 

W~(x) -> fc(x) 7 ,(l - l 5 )Mx)], (A.5) 

where Gy = Gp cosi?o 

The S matrix describing the transitions like p + p — > D + W + is defined by 

g = Te iJ d A x [C npD (x) + C npW (x) + C p c ^ pp {x) + . . .] ^ ^ A 6 j 

where T is the time-ordering operator and the ellipses denote the contribution of inter- 
actions irrelevant for the problem. 

For the evaluation of the effective Lagrangian £^ _>DW+ (x) we have to consider the 
third order term of the S matrix which reads 



S {3) = — J d 4 x 1 d 4 x 2 d 4 x 3 T([C npD (x 1 ) + £ npW (^i) + £ 



x [C QpD (x 2 ) + C npW (x 2 ) + C%+"(x 2 ) + ■■■} 

X [C npD (x 3 ) + £npw(x 3 ) + C^fa) + ...]) = 
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= -i I d 4 x x d 4 x 2 d 4 x 3 T(C^ PP (^)^nM^)C npW (^)) + ■■■ (A.7) 

The ellipses denote the terms which do not contribute to the transition p + p — > D + W + 
and the interaction £ npW (:r). The S matrix element S p p_,. DW+ describing the transition p 
+ p — > D + W + we determine as follows 

S pp-,dw+ = ~i J d' k x 1 d 4 x 2 d i x 3 T(£ p c ^ pp (x 1 )C npD (x 2 )C DpW (x 3 )). (A.8) 

For the derivation of the effective Lagrangian C Pl ^ DW+ (x) from the S matrix element 
Eq.(A.8) we should make all necessary contractions of the operators of the proton and 
the neutron fields. These contractions we denote by the brackets as 

(S p 3 pLdw + > = -*/ d 4 Xl d 4 x 2 d 4 x 3 (T(C^(x 1 )C a ^(x 2 )C npW (x 3 ))). (A.9) 

Now the effective Lagrangian £^~ +DW+ (x) related to the S matrix element (S p p_ >DW +) 
can be defined by 

(S p 3 p Idw + ) = i j 'd 4 x £ c P r DW+ (x) = 
= -if d 4 Xl d 4 x 2 d 4 x 3 (T(£ c P r PP (^i)A ip D(^)/:npw(a:3))). (A.10) 

In terms of the operators of the interacting fields the effective Lag rangian /3 P g (x) 
reads 

J d 4 xC PP T™ + {x) =-J d 4 Xl d 4 X 2 d 4 X 3 (T(£Pr PP (^l)^npD(^)A 1 pw(^3))) 

= -^C NN x (-igy) x (-g A ) Jd^x.d^Tiip^^pix^Dlix^W-ix,)) 
x (0\T(\p( Xl ) 7 a 7y(xi)][f(x 2 )7"n(i 2 ) - n c (x 2 )^p(x 2 )} [n(x 3 )Yl 5 p(x 3 )])\0) 
-ic NN x (-igy) x ( <?a) j d 4 x l d 4 x 2 d 4 x 3 T{[p c {x l )^p{x 1 )]Dl{x 2 )W;{x 3 )) 

x(0|T([p(x 1 ) 7 V(x 1 )]^(x 2 )7^(x 2 ) -nc(x 2 )Yp(x2)} [n(x 3 )Yl 5 p(x 3 )m- (A.ll) 

Since p + p — > D + W + is the Gamow-Teller transition, we have taken into account the 
W-boson coupled to the axial- vector nucleon current. 

Due to the relation n c (x 2 )^p(x 2 ) = —p c (x 2 )^n(x 2 ) the r.h.s. of Eq.(A.ll) can be 
reduced as follows 

J d 4 xC P r DW+ (x) = -J d 4 Xl d 4 x 2 d 4 x 3 (T(£ P r PP (^i)^npD(^)/:„pw(a:3))) 

= ^nn x (-igy) x g A J d 4 x 1 d 4 x 2 d 4 x 3 T(^(x 1 ) lal 5 p(x 1 )]Dl(x 2 )W-(x 3 )) 
x (0|T([p(x 1 ) 7 Vp c (^i)][p^^)7 M n(a; 2 )][n(x3)7"7 5 ^3)])|0) 
+C NN x (-igy) x g A J d 4 x 1 d 4 x 2 d 4 x 3 T(\£ c (x 1 ) 1 b p(x 1 )]D\(x 2 )W;(x 3 )) 

x (0|T([p(x 1 ) 7 V(x 1 )][p c (a; 2 )7' 1 n(x 2 )][n(a;3)7V^3)])|0). (A.12) 
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Making the necessary contractions we arrive at the expression 

J d 4 x £jg~* DW+ (x) = - J d 4 x 1 d 4 x 2 d 4 x 3 (T(C P ^~" pp (x 1 )C npB (x 2 )C npW (x 3 ))} = 

= 2 x C NN x (-igy) x g A J d 4 x x d 4 x 2 d 4 x 3 T([p c ( Xl ) 7a7 5 p(^i)] D U X ^ W~(x 3 )) 
x (-1) tr{ 7 V H)SM*i - x 2 )^(-i)S F (x 2 - x 3 ) 1 " 1 b {-t)S F {x 3 - Xl )} 
+2 x C NN x (-igy) x g A J d^^x^XsTilp-ix^pix^Dlix^W-ixs)) 

x (-1) tr{7 5 H)^(x! - x 2 )<f(-i)S F (x 2 - x 3 )Yl 5 (-i)S F (x 3 - Xl )}, (A.13) 

where the combinatorial factor 2 takes into account the fact that the protons are identical 
particles in the nucleon loop. 

In the momentum representation of the Green functions we get 

jd 4 xc*r DW+ (x) = 

— ^ n a Civm ——— [ d 4 x, f d4x ^ -ik, ■ (x, - x,) -ih* ■ (x 3 - x,) 

- i 9aC nn 8 ^ 2 J d Xl J ^ (2 ^ )4 e 

x T(\p°( Xl ) 7a7 5 p(^i)] Dl(x 2 ) W~(x 3 )) J a ^(k 2 , k 3 ; Q) 
_ io4 (- U 4 x, l d ^h.^^l e -ik 2 -(x 2 -x x )-ik 3 -(x 3 -x x ) 

igAC^^jdx^ (2?r)4 (2 ^ )4 e 

x T([^(xx) 7^1)] r>i(x 2 ) ^-(x 3 )) J""(A; 2 , A; 3 ; Q). (A.14) 

The structure functions J a ^ v {k 2) k 3 ; Q) and J^ v {k 2) k 3 ; Q) are defined by the momentum 
integrals 



fd 4 k ( a 



•' Y ^ - r f'f - 



rd 4 k r e 



M N - A; - Q + £; 2 M N — k — Q M N — k — Q — k 3 
J fM/ (k 2 ,k 3 ;Q) = 



J-^ : ^7" 1 ^— ^ ^1. (A.15) 

M N - A; - Q + A; 2 M N - fc - Q M N - A; - Q - A; 3 J 

We have introduced a 4-vector Q = ak 2 + bk 3 caused by an arbitrary shift of a virtual 
momentum with arbitrary parameters a and b. 

In order to obtain the effective Lagrangian £gg^ Dc+iyc (a;) of the transition p + p — > D 
+ e + + u c we have to replace the operator of the W-boson field by the operator of the 
leptonic weak current Eq.(A.5): 

Jd 4 xcir 1)c+Vc {x) = 

= i 9 K<hJl^ \d 4 x x [ d ^ d y e~ ik > ■ ^ - *0e-'*3 ■ (xs - *i) 

y v^87r 2 7 J (2tt) 4 (2tt) 4 

x T(^(xO 7a7 5 p(^i)] £>^ 2 ) foMxsMl - 7 5 )^c(a:3)]) J Q ^ 2 , A; 3 ; Q) 

39 



y v / 2 8vr 2 i V (2tt) 4 (2tt) 4 

x T{\p°( Xl ) 7 5 p(x 1 )] £)t (x 2 ) ftM* 3 )7,(l - 7 5 )^eOr 3 )D A; 3 ; Q). (A.16) 

Thus, the problem of the evaluation of the effective Lagrangian of the transition p + p 
— > D + e + + v c reduces itself to the problem of the evaluation of the structure functions 
Eq.(A.15). The momentum fc 3 is related to the 4-momentum of the leptonic pair. Due to 
the Gamow-Teller type of the transition p + p — > D + W + the contribution proportional 
to the 4-momentum of the leptonic pair turns out to be much smaller with respect to 
the contribution proportional to the 4-momentum of the deuteron k 2 . Therefore, without 
loss of generality we can set k 3 = in the integrand. This gives 

J aflu (k 2 ,k 3 ;Q) = 

d k I a 5 \ a \ v. 5 \ \ 



tr< 7" 7 ° s s ^7" s ^7*7 



vr 2 i {' ' M N — k — Q + k 2 ' M N - k - Q ' ' M N - k - Q . 

J^(k 2 ,h;Q) = 

^tr{ 7 5 , 1 s . 7" 1 -7V 1 4, (A.17) 

7T 2 z 1 M N -A;-Q + A; 2 M N - fc - Q M N - fc - Q J 

The evaluation of the momentum integrals at leading order in the 1/Mn expansion corre- 
sponding the leading order in the large Nc expansion due to the proportionality Mn ~ Nc 
in QCD with SU(N C ) gauge group at N c — > 00 [22] (see also Ref. [1]) yields 

J "" fa, h; Q) = 3 (kW - k^ a ) + + 2a) (fc^"" + fc^"°), 

^(fca, A; 3 ; Q) = <T4 M N J 2 (M N ) ~ 0(l/iV 2 ), (A.18) 

where the terms proportional to k 2 have been dropped, because they produce the con- 
tributions to the effective Lagrangian multiplied by d^D^x) vanishing by virtue of the 
constraint d^D^x) = 0. Then, J 2 (Mn) is a logarithmically divergent integral defined in 
the NNJL model in terms of the cut-off A D = 46.172 MeV such as A D < M N [1]: 



The cut-off Ad restricts 3-momenta of the virtual nucleon fluctuations forming the phys- 
ical deuteron [1]. Due to the uncertainty relation Ar A D ~ 1 the spatial region of virtual 
nucleon fluctuations forming the physical deuteron is defined by Ar ~ 4.274 fm. This 
agrees with the effective radius of the deuteron r D = 1/V^d^n = 4.319 fm. 

Keeping the terms of the same order in the large Nc expansion we get the structure 
functions 

J a ^(k 2 , k 3 ; Q) = 3 (Aftf - KgT) + 1(1 + 2a) (k^ + fc 2 T), 

J^(k 2 ,h;Q) = 0, (A.20) 
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The structure functions Eq.(A.20) define the effective Lagrangian £g^ De+ ^ c (a;): 



X 



(d.DKx) - d u Dl(x)) + 1 (1 + 2a) (d.DKx) + d v D\{x))\ . (A.21) 

In order to remove uncertainty caused by the shift of virtual momenta of the one-nucleon 
loop diagrams we demand the invariance of the effective Lagrangian £^^ Dc Uc (x) under 
gauge transformations of the deuteron field 

Dl(x)^Dl(x) + d,n(x) : (A.22) 

where Q(x) is a gauge function. The requirement of the invariance of the effective La- 
grangian of the low-energy transition p + p^D + e + + z/ c under gauge transforma- 
tion Eq.(A.22) imposes the constraint a = —1/2. This reduces the effective Lagrangian 
£efT De+l/e (z) to the form 

£ p r De+, o(x) = g A C^^ Dl{x) ip<(x)Yl 5 p(x)} $*(x)>f(l ~ T^eW]- (A.23) 

This effective Lagrangian is defined by the structure function J^^iki-, k 2 ] Q)'- 

J a » u (k 2 , h; Q) = 3 {kW - k^n- (A.24) 

The structure function J" a ^(/c 2 , k 2 ; Q) does not depend on the mass of virtual nucleons 
and according to Gertsein and Jackiw [37] can be valued as the anomaly of the one-nucleon 
triangle AAl^-diagram. The requirement of gauge invariance applied to remove ambigu- 
ities of the structure function J a ^ 1 '{k 2 , k^, Q) and to fix the contribution of the anomaly 
of the one-nucleon loop AAV-diagrams is in complete agreement with the derivation 
of the Adler-Bell- Jackiw axial-vector anomaly performed in terms of one-nucleon loop 
Al^F-diagrams (see Jackiw [54]). 

The effective Lagrangian £^ D_>c+nn (a;) describing the low-energy transition P e + D — > 
e + + n + n can be obtained by the way analogous to Cl^ Dc+Uc (x) and reads 

£ pjWn„ (x) = _^ A c NN ^^^(x) [n(x)n 5 n c (x)] foUsMl -T 5 )^*)]- 

(A.25) 

In the low-energy limit when (x) [p c (x)'~f^'~f 5 p(x)] — > — 2i M^Dl(x)\p c (x) , y 5 p(x)} the 
effective Lagrangian Eq.(A.23) reduces itself to the form 

£ir Be+Ue (x) = -i 9 AM N C NN ^^Di(x) mx)^p(x)} [^ c (a;)7^(l -7 5 )^e(a:)]. 

(A.26) 

The low-energy reduction D fJil/ [n(x)^^ 5 n c (x)] — > — 2i M^D u (x)n(x)^n c (x) applied to 
the effective Lagrangian nn (x) gives 

£^e + nn (x)= .^ MnCnn ^ (A.27) 

41 



The effective Lagrangians Eq.(A.23) and Eq.(A.25) as well as Eq.(A.26) and Eq.(A.27) 
testify distinctly that the transitions p + p^D + e + + z/ e and u e + D — ► e + n + n are 
governed by the same dynamics of low-energy nuclear forces in agreement with charge 
independence of weak interaction strength. 

For the evaluation of the effective Lagrangian ~* Ve1xp (x) of the transition v e + D 
— > v c + n + p in the NNJL model one should use the following Lagrangians 

4 p d(^) = -igv\p(x)^n c (x) -n(x)^p c (x)]D^x), 

C^ np (x) = C NN {[p(x) Yl 5 n c (x)] [n^x) J^J 5 p^)} + W) 7V(x)] [n^x) ^p(x)]}, 
£nnzO) = gx \p(x)j u j 5 p(x) - n(x)^ u ^n(x)} Z u (x), 

Z v {x) || i^ c (x)Y(l - 7 5 )^e(x)]. (A.28) 

The effective Lagrangian C^ > ~* Uenp (x) of the low-energy transition u e + D — > v c + n + 
p is then defined by 



V2 07T 



(x) = -SaCm^^(i) \p(x)n b n c (x)} [^(^/(l-T 5 )^^)]. 

(A.29) 

In the low-energy limit when D fiu (x)\p(x)'y fl '~y 5 n c (x)} — > —2iM^D u (x) \p{x)rfn c {x)] the 
effective Lagrangian reduces itself to the from 

C^W^jAM^NNl^^W^l^WHif^ll-flif,)]. (A.30) 

For the evaluation of the matrix element of the transition z/ c + D^z/ e + n + p 
one would use the wave function of the np pair in the standard form \n{jpi)p{jp 2 )) = 
at(pi, ffi) a>p(p2, cr 2)|0). The contribution of low-energy nuclear forces to the relative 
movement of the np pair in the 1 S state should be described by the infinite series of 
one-nucleon bubbles evaluated at leading order in the large Nc expansion. The result 
should be expressed in terms of the S-wave scattering length a np and the effective range 
r np of low-energy elastic np scattering in the 1 S state. 

Now let us obtain the contribution of the nucleon tensor current Eq. ([4.10| ). In terms 
of the structure functions the effective Lag rangian 5CJ^ c (-^) 



j d 4 x6C^ De+Uc 



x) 



CmGvJt_ f ,4 f d*x 2 d*k 2 d 4 x 3 d 4 h -ik 2 ■ (x? - x^-ifa ■ (x 3 - xi ) 
^ A 87r 2 v / 2 2M N 7 V (2tt) 4 (2tt) 4 

x T(\p^ Xl ) 7a7 5 p(^i)] DUx 2 ) [^ c (x 3 ) 7 a(1 - 1 5 )M^)}) J^Ck, k 3 ; Q) 



9 a 



^ r d A x 2 d A k 2 d A x 3 d% -ik 2 ■ (x 2 - Xl ) e -ik 3 ■ {X 3 - Xi) 



87r 2 v^2M N y l J (2tt) 4 (2tt) / 
x T{\p(xi) 7 5 p(x 1 )] Dl(x 2 ) fe e (x 3 )7A(l - 7 5 )^e(x 3 )]) J^\k 2 , k 3 ; Q). (A.31) 
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The structure functions are given by 

f „ K 1 1 \ pr 1 



= / — tr<J 7 a 7 5 — ? A-^^TT f a7V 



n z i { M N — k — Q + k 2 M N - fc - Q M N - fc - Q - k 3 

J^\k 2 ,k 3 ;Q) = 



fd 4 k( 5 



' ^(7^ 1 ^ 7 V ^— ^ ^1, (A.32) 



M N -k-Q + k 2 M^-k-Q M N - k - Q - k 3 

where a 4- vector Q = a k 2 + b k% is an arbitrary shift of a virtual momentum with arbitrary 
parameters a and b. 

The evaluation of the structure functions Eq.(A.32) at leading order in the large Nq 
expansion gives the following effective Lagrangian 5C V J^ ° Ua (x) caused by the contribu- 
tion of the nucleon tensor current 

5£ c P r Dc+ -(x) = g A C NN {Dl(x) [f(x) 7 " 7 yx)] - ^- PD^x)^)^)]} 

x[^ e (x) 7 l/ (l-7 5 )^e(x)]. (A.33) 

It is seen that the coupling constants of the effective Lagrangian depend on the arbitrary 
parameter a caused by a shift of a virtual momentum. 

The analogous expression one can get for the effective Lagrangian of the transition V e 
+ D^e + + n + n caused by the contribution of the nucleon tensor current as well 

5£lf-° + ™(x) = g A C NN ^ {2V(s) [n(x)>f>fn°(z)] - J^L d»D„„(x)[n(x)^n c (x)}} 

x[^ e (x)7 I/ (l-7 5 )^e(^)]. (A.34) 

In the low-energy limit when 

DU X ) [P c (xh^ 5 p(x)] ^-2tM^Dl(x) [pc(xh 5 p(x)}, 

D^(x) [n{x)Yl 5 n c {x)] -> - 2 1 M N D u {x) [n(x) 7 5 n(x)] (A.35) 
the effective Lagrangians Eq.(A.33) and Eq.(A.34) can be recast into the form 

Cr c+m W = - l gAM^C^^^iD u {x)[n{x) 1 "n c {x)] [^(x)-f(l - 7>e(*)], 

(A.36) 

where £ is an arbitrary parameter related to the parameter a as follows 

£=1(4 -a). (A.37) 
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The total effective Lagrangians of the transitions p + p^D + e + + z/ e , z/ e + D^e + + 
n + n and z/ c + D^z/ c + n + p are defined by 

C:f^(x) = i (1 + 9a M N C NN D v {x) [p(xh 5 n c (x)] ~ 7 5 )<M*)], 



(A.38) 

where £ is obtained by using the relation g T = ^3/8 gy and is defined by 

e=\/|e (A.39) 

Under the assumption of isotopical invariance of low-energy nulcear forces the best agree- 
ment with the recommended value for the astrophysical factor for the solar proton burning 
[29] and the contemporary experimental data [36] on the cross sections for the anti- 
neutrino disintegration of the deuteron z/ c + D^e + + n + n and z/ c + D^z/ c + n + p 
caused by charged and neutral weak current, respectively, we obtain at £ = 0. 
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Figure caption 



• Fig.l. One-nucleon loop diagrams of the contribution of the effective coupling 
[p(x)'y 5 n c (x)}[n c (x)'y 5 p(x)] to the effective Lagrangian of the Ml transition n + p — > 
D + 7. 

• Fig. 2. One-nucleon loop diagrams of the contribution of the effective coupling 
[p(x)^ a ^n c (x)][n c (x)^ a ^p{x)] to the effective Lagrangian of the Ml transition n + 
p -> D + 7. 

• Fig. 3. One-nucleon loop diagrams describing the effective Lagrangian C^* Dc+Uc (x) 
of the low-energy transition p + p^D + e + + z/ e . 
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